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ABSTRACT 

The radiative- energy current due to line radiation is calculated 
235 

in a U plasma over a temperature range of 5000“K to 800G®K. Also 

12 2 

a variation in the neutron flux of 2x10 neutrons/ (cm -sec) to 
16 2 

2x10 neutrons/ (cm -sec) is considered. The plasma forms a cylinder 
with a diameter and height of one meter. 

To calculate the radiative- energy current, a knowledge of the 
atomic state densities within the plasma is necessary. A rate equation 
formalism is developed to solve for the atomic state densities along 
with a model for the energy levels in neutral and singly-ionized 
uranium. 

Because the electron states in uranium lie below 5eV, recombination 

is the principle excitation mechanism. At and above 6000*K, inversions 

were found and at all temperatures the line radiation at line center 

was greater than the corresponding black-body radiation. An example 

of this is the 28763 cm~^ to 5762 cm”^ transition in neutral uranium, 

— 6 2 

where the Planck function at 5000®K is 6.49xl0~ ergs/cm and the 

-4 2 

calculated radiative-energy current was 1.49x10 ergs/cm . There are 

negligible differences in the radiative -energy current at 6000" K for 

variations in the neutron flux. Also the average opacity, which varied 
5 2 

from 100 to 10 gm/cm , over the frequency range of line radiation is 


calculated. 
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CHAPTER I. DISCUSSION OF THE PROBLEM 

In a gaseous core reactor non-equilibrium radiation will be emitted 
from the plasma due to the presence of a non-Maxwe Ilian electron flux and 
rapid recombination caused by the relatively low electron temperature. 

The frequency distribution and intensity of the emitted radiation are 
important because of the potential applications for energy conversion and 
chemical processing as well as the material constraints they put on the 
reactor design. This study is concerned with the line radiation emitted 
over a range of operating conditions envisioned for a gaseous core uranium 
reactor. ^ ' Thus a uranium plasma at a pressure of one atmosphere with 
temperatures ranging from SOOO®K to 8000«K, (the boiling point of uranium 
being 4407®K at 1 atm) will be investigated. 

Normal treatments of the problem are insufficient because of the 
physical parameters of this unique plasma. A "coronal" model, where the 
principal method of depopulating excited states is through spontaneous 
decay, is inadequate because of the large collisional rates due to the 
large electron population. An equilibrium model, where the collisional 
processes completely dominate state relaxation, is also inappropriate 
because the large non-equilibrium excited state densities lead to non- 
negligible radiative relaxation from excited states. Hence, the state 
densities are explicitly calculated using a rate equation approach 
and then the emitted radiation is evaluated. 

The plasma conditions cited above are characteristic of a subcritical 
uranium plasma. Critical reactors have higher tenq)eratures (center line 
temperatures of 40,000®K) and pressures (approaching 500 atmospheres ^^^ ) . 
However, the results of these calculations should be applicable to the 
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outer layers of a critical reactor. Also, this temperature range encom- 
passes the 6000*’K temperature of the proposed 5 Megawatt Los Alamos 
experimental reactor. 

In addition to the effect of a temperature variation, a change of 
the neutron flux of 2 x 10^^ to 2 x 10^^ neutrons/ (cm^-sec) is considered. 
Finally, the effect of the size of the plasma on the emitted radiation 
is studied. The sizes of interest range from a cylinder 4-cm diameter 
and 18-cra long^^^ to a larger size corresponding to the Los Alamos experi- 
mental plasma of 105.5-cm long and 102.4 cm in diameter. 

A. Previous Work Related to the Problem 

NASA has been interested in a gaseous core reactor for space pro- 
pulsion purposes for approximately ten years. Much of the early work 
was carried out by workers at United Aircraft who assumed the plasma 
emission had a black-body spectrum. The first detailed calculation for 
uraniiiin was carried out by Parks et al ^^^ at General Atomic. The under- 
lying assumption in their work was that equilibrium statistics can be 
used to determine state densities within the plasma. However, recent 
studies have shown that a gas in the presence of a neutron-induced 
source will not have a Boltzmann distribution of excited states. This is 
because a non-Maxwell ian electron flux is produced within the plasma. Also, 
at the high ion and electron densities, recombinational excitation becomes 

. ri2‘) 

an important mechanism.'- ^ " 

An early treatment of radiation transport in a plasma was done by 
Holstein in which he considered trapping of resonance radiation. 

His results are significant because he showed the importance of considering 
the frequency dependence of the absorbing medium. However, an exact 
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calculation of the density of excited states and non-equilibriian excitation 

mechanisms were ignored In his analysis. A later analysis assumed an 

excitation temperature for a transition and then used a thermodynamic source 

term for the emitted radiation. Poraraning^^^^ formulated the source term 

through a rate equation formulation, and then did some sample calculations 

in a slab geometry. Later the concept of multi-temperature electrons was 
fl7 18T 

used by Eddy^ ■ ^ where a finite number of excitation temperatures was used. 

In the present analysis a rate equation approach is used to calculate 
the excited state densities. Unlike previous calculations, the radiative- 
energy intensity is evaluated and then substituted into the rate equations. 
This allows state densities to be calculated when neither collisional nor 
radiative processes dominate relaxation within the plasma. Also, state 
densities can be calculated when either collisional or radiative processes 
dominate relaxation. From these state densities the radiative -energy 
current is calculated. Another feature of the present analysis is that it 
can treat a non-uniform plasma (i.e., variations in temperature, pressure 
or flux) . 

B. Physical Applications 

Well referenced discussions of the uses of gaseous core reactors 

fl9 20) 

exist in the literature, ^ ^ hence an extensive list of references is 

not included here. 

The initial impetus for research on a gaseous core reactor was to 
develop a nuclear rocket engine. Shorter transit time for space flights, 

(80 days for a round trip to Mars as opposed to 200 days using conventional 
propulsion systems), and savings in weight, (up to an order of magnitude 
over chemical propellants on long flights )j, make nuclear .rocket engines 




desirable. Another advantage of the nuclear rocket is that the propulsion 

system can be adapted to simultaneously provide power for the space craft 

after take off. In this application the emitted radiation is important 

because it is used to heat the propellant. 

Several characteristics of the gaseous core reactor also make it 

attractive for terrestrial power production. Thermodynamic cycle 

efficiencies of approximately 65% can be attained in a plasma core-MHD 
( 21 ) 

device at 5000® K. Because of the lower densities in a gas core 

reactor, smaller inventories of nuclear fuel would have to be maintained 

at the plant site thus reducing production costs. This reduced inventory 

would also decrease the amount of radioactive waste at the plant site. 

Finally, because the fuel is a gas, refueling could be done continuously, 

eliminating the need for lengthy reactor shutdowns. 

In addition to normal power plant uses, because of its high thermal- 

neutron flux-density, the gas core reactor could be used as a breeder of 

U from thorium, Non-therraal energy- conversion cycles could be developed 

utilizing the non-Maxwell ian electron flux and fast recombination rate for 

(22 231 

direct nuclear pumping of lasers. * Further, radiation at a desired 
frequency could be produced by the selection of an appropriate seed gas 
to be used in chemical processing. 

Any of the processes dealing with light, lasers, or chemical processing 
require a detailed knowledge of the emitted radiation energy spectriun and 
intensity. Also in any gaseous core configuration the radiation mast be 
transmitted through a semitransparent wall. Absorption in this "window*', 
which 0 an cause extensive heating, is frequency dependent. These consider- 
ations motivated the present study. 
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C. Terminology 

Several terms, namely gain, directed radiative-energy intensity, 
total radiative-energy intensity and radiative-energy current are defined 
here. 

The gain, y, describes the amplification, (or absorption), of photons 
in the plasma. This can be expressed 






8eS 

e t 

3Cmc^)lcT 



where is the spontaneous transition probability, c, the speed of 
light, V, the frequency of the transition, g(v), the line shape, N^, the 
upper state density, N., the lovjer state density, N , the density of 
electrons within the plasma, N^, the density of ions within the plasma, 
g j^/gj » the ratio of the upper state to lower state degeneracies, e, the 
electron charge, m, the electron mass, k, Boltzmann’s constant, h, 
Planck's constant, and T is the temperature of the plasma. The units 
of gain are cm ^ conforming to the CGS units used throughout this work. 

If the plasma is homogeneous, the intensity at a distance s from 
the origin due to an attenuated intensity at the origin can be written as 

I(s,v) = I(0,v)exp(Y(v)s) . (1-2) 

The directed radiative- energy intensity, I.^(e,(j)) is the energy per 
unit time per unit area per steradian per unit frequency at frequency v. 
From this qv,antity the total radiative- energy intensity, F^(s), can be 
calculated by integrating the directed radiative-energy intensity over 
all angles as follows: 
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F^(s) - f i^(e,<)))dn. (1-33 

where the units for F^(s) are the same as 1^(0, (|>3 With the angular de- 
pendence integrated out. The total radiative- energy intensity combined 
with the Einstein B coefficient for absorption or stimulated emission 
gives the reaction rates for these processes within the plasma and 
can be written 

R(s) = I BF^(s). (1-4) 

To calculate the radiative energy- spectrum emitted at the surface 
of the plasma, the radiative-energy current, J^, must be known. This is 

^ At 

obtained by multiplying the directed radiative energy flux times 

At 

where is the direction the current is being evaluated and integrating 

At 

this quantity over all Hence 

= f I^(6,<|))r*ft dQ' . (j_5) 

This result has the same dimensions as the total radiative-energy 
intensity. 

In Chapter II a theoretical formalism for rate equations that describes 
the state densities within the plasma is discussed along with the radiation 
transport problem. The gain of Eq. (I-l) and the total radiative-energy 
intensity of Eq. (1-3) are used in Chapter II. Chapter III describes the 
computer code developed using this formalism and Chapter IV describes the 
plasma and the physical quantities that are used in the analysis. In 
Chapter V the results of these calculations are presented and discussed. 
Chapter VI contains the conclusions drawn from this work and suggestions 


future research. 



CHAPTER II. RADIATION TRANSPORT PROBLEM 


A. Physical Parameters 

Electron-ion recombination and non-Maxwell ian electron excitation 
lead to non-equilibrium radiation emission by a plasma. This study 
deals with the radiation emitted by such a plasma. The plasma processes 
that lead to the emission are illustrated in Fig. (II-l). The electron 
flux is produced by fission fragments within the gas. This flux in turn 
excites the constituents of the plasma leading to line radiation. The 
radiation is absorbed or amplified by the medium leading to a shift in 
the excited-state densities of the gas. This shift in excitated state 
densities will be significant when the line radiation they produce is 
large enough to contribute to the rate equations. The interplay between 
the radiation and state dev.sities within the plasma results in a non- 
linear problem. 

To solve for the emitted radiative current, the excited-state den- 
sities within the plasma are first calculated in the absence of any radia- 
tion and then the radiation intensity is computed. The state densities 
are recalculated using the computed radiative intensity. Then the 
radiative intensity is evaluated with the new state densities. Tliis 
process is repeated until the solution for the radiative intensity con- 
verges. In this manner the nonlinear problem is solved. 

In this work a cylindrical plasma is assumed because this most 
nearly models the configuration planned for the Los Alamos gaseous-core 
experiments as well as the design proposed by United Aircraft for 
prototype models of gaseous-core rocket engines. In this analysis the 
plasma is broken up into concentric rings. This is done so any temperature. 
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Figure II-l. Flow diagram of physical processes within the 
plasma. 
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flux or density gradients within the plasma can be included. Within a 
given region all physical properties of the plasma are assumed to be 
constant but these properties change from ring to ring to account for 
gradients. This model presupposes constant physical properties along 
the axis of the cylinder. 

B. Rate Equations 

3 ’) 

The rate equations are a system of equations * that describe the 
atomic-state densities within the plasma. They have the form 


E C . 

pq p q 


(II-l) 


where Np is the state density of the p^^ state of the atom, is the 

collisional-radiative matrix coefficient and is a collisional constant 
describing recombination. Figure gives a schematic representation 

of the processes that determine the coefficients in the rate equations. 

C denotes collisional processes between states, i.e. excitation and de- 
excitation rates. . R denotes radiative processes. Both radiative and 
collisional processes will take electrons from the upper to lower state 
and from the lower to the upper state. Hence, the paths in Fig. (II-2) 
have arrows at both ends but each path is labeled only once for simplicity. 
In this work the first subscript of the matrix element is the designation 
of the initial atomic state and the second subscript is that of the final 

atomic state. 

Photon absorption is the process that takes electrons from a lower 
state to an upper state while spontaneous and stimulated emission account 
for the radiative processes that take electrons from an upper state to a 
lower state. Stimulated emission is when a photon with energy equal to 
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the difference in the state energies collides with an excited atom 
causing it to de-excite emitting a photon in phase with the original 
photon. Ionization causes bound electrons to become part of the continuum 
while recombination is the inverse process. Since recombination takes 
free electrons from the continuum into atomic states in the plasma it is 
represented by the right hand side of Eq. (II-l) . The states used in this 
work are described in the model in Chapter IV. In uranium three-body 
processes dominate ionization and recombination thus no radiative channels 
are seen from any bound state to the continuum in Fig. (II-2) . Also, as 
will be described in Chapter IV, there are no radiative processes within 
the families of upper states and lower states as illustrated by no 
radiative processes between the ground state and state p. 

Processes that are considered in the present study are listed below 
and discussed in Chapter IV. 

1. Spontaneous Emission of Radiation: 

X(q) XCp) + hv, (n-2) 

where p and q represent different energy levels and V is the fre- 
quency of this transition. The spontaneous transition probability 
is A with dimensions of sec“^, and is equal to the dipole moment 

qp 

for the transition. 

2. Absorption of Radiation: 

X(p) + hv XCq), Cn-3) 

The rate for this interaction is p(v) where p(v) is the radiation 
density at v per unit V and is the Einstein B coefficient. V/here 

B has dimensions of cm^~sec ^-ergs and p has dimensions of ergs-sec- 

-3 
cm . 
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3. Stimulated Emission of Radiation: 
X(q) + hv -► XCp) + 2 hv. 

The rate for this interaction is p(v) B 


(I 1-4) 


iff* 


4. Electron Excitation: 

X(p) + e XCq) + e. (II-5) 

The rate for this process is given by Ja^pCE) 4>g(E) dE where CTpq(E) 

is the excitation cross section and <|)-(E) is the electron flux per 

6 

- 2-1 

unit energy. The flux has units of cm -sec and the cross sections 
has units of cm . All cross-sections are zero for energies below 
threshold. These cross sections are for all states including the 
ground state. 

5. Electron De-Excitation: 

X(q) + e X(p) + e . (II-6) 

The rate for this process is given by / o _(E) <j>_(E) dE. 

qp ® 

6. Ionization of State q: 

X(q) + e -*■ X* + 2e. (H-7) 

The rate for this process is given by <l>gCE) dE where Oqj^(E) 

is the ionization cross- sect ion. 

7. Recombination into State q: 

i) Two-body recombination 

X"*" + e X(q) + hv. (H-8) 

ii) Three-body recombination 

+ 2e -*■ X(q) + e . (H-9) 

The total Ci*e., combined) rate for these processes can be writteh 
jn^CE) lij^CE) ot^ff^^^ dE where n^CE) and n^CE) are the electron and ion 
number densities at E per unit E and ®^g£fCE) is the "effective** two-body 
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(2 3") 

recombinational branching coefficient to state q'- ’ The effective 
two-body recombinational branching can be written as 

a^ff(E) = a^CE) + n^(E) a^CE), (II-IO) 

where is the two-body recombinational branching rate to state q and 
is the three-body recombinational branching rate to state q. The branching 
ratio is the fractional part of the recombination coefficient that accounts 
for recombination into state q. 

All of the above processes are taken into account for the constituents 
of the uranium plasma. If another plasma is to be considered, these pro- 
cess will have to be re-examined to see which ones are important. In this 
analysis, atom-atom and atom- ion collisions have been neglected because 
they do not contribute significantly to excitation and ionization at thermal 
energies due to the small cross sections in this energy range. Direct 
fission fragment excitation is also neglected due to the low fragment 
density. 

When the processes described by Eq. (II-2) through Eq. (11-9) are 
combined and an infinite medium is assumed so diffusion losses can be 
neglected, the rate equation obtained for state q becomes: 

Z V % * E “pq’ "p 

p<q Fq 

- V ”q * E V "p 

p>q P<q 

vM p?q 

- N^ya^iCE) 4>g(E)dE + y* n^CE) n^^ (E) a2£f(E)dE. (11-11) 
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To describe specific plasma conditions , the sum over all state 
densities of a given atomic species should equal the total number of atoms 
of that species. This can be written: 




ClI-12) 


The total radiative-energy intensity, F^(0), enters into the rate 
equations because it is proportional to p(v). This proportionality is 


p(v) = 


F^(0) 


(11-13) 


here F(v) has units of ergs/cra"^. 

Equation (11-11) can be simplifed in the case of a steady state. 
When 4^ is set equal to zero, Eq. (II-ll) has the form of Eq. (II-l) 


where 


\ "i(E) a2f^(E) dE, 


(11-14) 


(Apq + P(V) Bp^) ^y^pqCE) <^.^(E)dE, 

■ £ % * P'"’ V -£/qp' 

p<q . 


qp(E) *JE)dE 


-/c 


<(r^CE)dE - 




PCv) b_. 


p > q 


p = q 




With the coefficients in Eq. (11-15), Eq. (II-l) can be solved to give 
the atomic state densities of the various energy levels in the plasma 
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under consideration. Once these state densities are known, the emitted 
radiation from the plasma is calculated. 

C. Radiation Transport 

1 . Plasma Source and Attenuation 

The volumetric source of radiation per unit frequency per unit 
time within the plasma is evaluated utilizing the atomic state densities 
calculated from the rate equations. This source has two components and 
can be written: 


f(v) = fj^Cv) + fp(v). 


(11-16) 


where f. is the source of line radiation due to atomic transitions and 
L 

f is the source of continuum radiation due to bremsstrahlung within the 
P 

f24‘) 

plasma. The first terra is given by'- 


AN hv g(v) 


CII-17) 


where A is the spontaneous transition probability, h is Planck's constant, 
V is the frequency of the transition, g(v) is the line shape and is 
the atomic state density of the upper level of the transition calculated 
from the rate equations and 6(V-V ) is the delta function. Units for 


the volumetric source are ergs/cm . 

The continuum source is given by 
6 


(25) 


f = 


4eS N. / 2 /, \ 

4/inc__ exp J 

.2 4 V 2-jrhv \ kT/ ’ 

3TTm c ’ 


(11-18) 


where e is the electronic charge, m is the mass of the electron, c is 
the speed of light, k is Boltzmann's constant, T is the temperature of the 
plasma, N^ is the electron density in the plasma, and is the ion density 
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of the plasma. Equation (11-18) neglects contributions due to free-bound 
transitions (recombination events). The photons emitted during recombina- 
tion events represent a small contribution to the source at the frequencies 
considered in this work and this allows the plasma emission to be of the 
form of Eq. (11-18).^^'^^^ 

The attenuation at a given frequency per unit frequency per unit 
length for this radiation can be written 


Y(v) = + Yp(v) , 


(11-19) 


where Yt is the attenuation due to the atomic-state densities within 
L 


the plasma and Yp is the plasma absorption coefficient. 


(25) 


The atomic 

attenuation is due to bound-bound transitions in the atoms while the 
plasma absorption coefficient is due to free-free transitions. Free- 
bound transitions are neglected because their contribution to the attenua- 
tion are negligible. The attenuation due to atomic-state densities can be 
written as: 







«i 

— N. 

gj 3 


(11-20) 


where is the spontaneous transition rate (Einstein A coefficient), c 

is the speed of light, g(v) is the line shape, is the state density of 

the upper state, Nj is the state density of the lower state and g^^/gj is 

the ratio of the degeneracies of state i to state j. The corresponding 

(25) . 

expression for the plasma absorption ^ is 


8eS N. 
e 1 

3(mc^)kT 



(11-21) 
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Ones the source rs.te snd 3.ttenu3.tion of xedietion sre known within & 
the plasma, the directed radiative-energy intensity along path length s 
can be calculated. 

2 General Relationships for the Total Radiative-Energy Flux and 
the Radiative-Energy Current 

To do these calculations an expression for the directed radiative- 
energy intensity must first be derived. Figure (II-3) shows a cylindri- 
cal volume element, along path length s, that has surface area dA normal 
to s and length ds. Within this volume element there is a volumetric 
source f(v,s), corresponding to the source in Eq. (11-16) and an attenu- 
ation of radiation, Y(v,s), through the element. The divergence of the 
direct radiative -energy intensity is defined by dJ2 about Jl, where is 
the angle s makes with the origin. 

The energy incident upon the surface can be written as 

I(v,s) dV dt dA, (11-22) 

where I(v,s) is a directed radiative- energy intensity, dv is the frequency 
interval, dt is the time interval and s = s((2). Similarly the radiative- 
energy intensity leaving the volume is I(v,s) + dl(v,s) and the energy 
leaving is 

tI(V,s) + dI(V,s)] dv dt dA. (11-23) 

Inside the volume element, the radiative- energy intensity will be 
incremented due to the volumetric source, which can be written as 

£(V,s) dV dv dt, (11-24) 

where dV = dA ds. The radiative- energy intensity will be attenuated 
along path length s, which can be written as 

(11-25) 


I(v,s) y(V,s) ds dv dt dA. 



Figure II-3. Cylindrical volume element with surface area dA and 
length ds. 
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Combining Eqs. (11-22) through (11-25), the expression for the intensity 


becomes 


Y(v,s) I(v,s) + £(v,s). 


(11-26) 


Equation (11-26) can be integrated over the path length s to obtain the 
integral form of the transport equation 


I(v,s) = I(v,o) expj / Y(v,s')ds' 

\ ) r ^ 

+ exp / Y(v,s')ds' / f(v,s') ds’ exp -I Y(v,x*)dx' 


Equivalently Eq. (11-27) can be rewritten to express the directed radiative- 
energy intensity at the origin when the directed radiative-energy 


intensity at s is zero as 


,s')ds» 


I(v,o) = exp /Y(v,s*)ds’| 

I 

J' f(v,s') ds' exp I ^(v. 


x) dx 


(11-28) 


In Appendix A expressions for the radiative-energy current and the total 
radiative- energy intensity will be derived in terms of the directed 

radiative- energy intensity. Analytic expressions for these quantities 

will be obtained for two simple geometries. The general expressions 

from Appendix A can be expressed for a cylindrical geometry as follows: 

When 0 is the azimuthal angle and ^ is the angle made with the x-axis, 

the radiative-energy current at the origin in a direction is written 

- 2tt .1 


/• 

>,fi») = J J 


d(cos6) (a*fi')I^(0,<{)), 


(11-29) 
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where and fl* are unit vectors in their respective directions and 


I(v,s[e,(|)l). (11-30) 

Also the total radiative-energy intensity can be written 


F^(o) » 


^2tt a 

ii: 


d(cose) 1^(0, 4>), 


(11-31) 


These expressions will be used in the subsequent developnent of the 
radiation transport problem. 


3. Radiation Transport in a Cylinder 

Since the total radiative-energy intensity is used in the rate 
equations for each region* Eq. ill-31) will form the basis for the 
following analysis. 

In order to handle any density and temperature gradients involved, 
the cylindrical plasma is divided into regions within which the state 
densities and volumetric source are approximately constant (Fig. II-4). 
The simplest case is when the plasma is homogeneous (i.e., there are no 
gradients in the electron flux or density of the plasma and it can be 
characterized by a single temperature) . In this case the plasma can be 
described by two regions. The inner region, whose radius is half the 
plasma radius, so the average state densities within the gas can be 
calculated, will account for the radiative- energy intensity to be used in 
the rate equations. The calculation employs an iterative technique. 

The convergence of the total radiative-energy intensity at the surface 


implies the convergence of the radiaitve-energy current, which will ulti- 
mately be calculated. 

If the plasma is non-homogeneous, then the number of regions is 





Figure II-4. 


Cross-sectional cut of a cylindrical plasma that is divided 
into four different regions. Excited state densities are 
acciimfir) to hft constant over each region. 
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determined by the gradients. For example, if the core region is at 
7000"K and the outer region is at 6000* K, then there would be three 
regions. The inner region would describe the core while the outer two 
regions would allow the calculation of the state densities and emitted 
radiation for the 6000*K region. This process of subdividing the plasma 
can be repeated to obtain a given accuracy for larger gradients in 
physical properties, as in the cases of the actual rocket engine config- 
uration where the ten 5 )erature goes from a value of 10Q,000*K at the 

(19) 

center to 10,000*K at the plasma surface. 

To calculate the total radiative -energy intensity emitted from the 
plasma an appropriate boundary condition on the outer surface is required. 
Since only the emitted radiation is of interest, the appropriate initial 
condition on I(v,s) is that the inward directed radiative-energy flux is 
zero. This condition can be written as 

where s indicates the path originates on the surface of the plasma. In 
max 

the method used in these calculations, the origin is shifted from region 
to region to evaluate the total radiative-energy intensity within each 
region. Hence, in the following analysis F^(0) is of central importance. 


Now the geometry of the problem must be considered. First the 
boundaries of the cylindrical shells must be known (Fig. II-5). The equa- 
tion for a cylindrical shell of radius r^ with its central axis displaced 
along the x-axis by a distance r^ is 


2 2 

T sinQ - 2 


2 2 

r r.sinG cos(|> + (r -r.) » 0. 

i ^ 


(11-33) 



Figure II-5. Cross section of cylinder with origin of coordinate system 
displaced by a distance -r^ along the x-axis. Here the 
z-axis is perpendicular to the plane of the slice shown. 
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In Fig. (II-5) A = 180® - (j) and 6 is not shown because the z-axis is 
perpendicular to the plane of the slice shown. Equation (1I«34) can be 
solved for the roots of r to give 


r.coscb + Vr? - r?si 

r = -i - V J 1 

+ sin0 


2 

Sin (|i 


(11-34) 


where r^ is the longer distance and r_ is the shorter distance to the 
cylindrical shell from the origin for a given 6 and Another quantity 


of important is the angle (Fig- II-5). 


a s 


or, solving for A 


13 


sin A. . = 

ij r. » 


Aij = (77) 


Sin A^j can be written 


(II- 35a) 


(II-35b) 


Because the mathematical form of the path length varies when the 
path tranversed by a ray originates on the top and the side of the cylinder, 
two explicit forms of 1^(6, <|)) must be obtained. In this formulation of 
the problem all calculations will be made for the plane that cuts the 
cylinder in half in the x-y plane. The length of the cylinder is 2b and 
its radius is a (Fig. II-6). To solve the problem, two different volumes 
must be considered. The forward volume to the right of the origin, 

[ IT TTl 

- -j , — ^ and the backward volume to the left of the origin, 
where (|> is the angle measured from the x axis in the x-y plane. This 
division is made because in the backward volume a path s can cross the same 
cylindrical boundary twice while in the forward volume the crossing will 
only happen once (Fig. II-7). With this division of the cylinder, Eq. 
(11-31) can be written 


O 9 

;ln. 
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F^(o) = F^(o) + F’(o) (11-36) 

F*(o) and F~(o) are the forward and backward total radiative -energy 
intensities defined by 

7T/2 1 

F^(o) = /* d(j> /* d(cose) I^(e,(t)), (II-37a) 

- ff/2 *'-1 

and 

3tt/2 1 

F“(o) = f d<l> f d(cos0) I (e,<|)) . (II-37b) 

Tf /2 *^1 

The + or - sign indicate the integration is done over the forward or 
backward volume of the plasma. Let the plasma be divided into n regions 
and let F.^(0) be evaluated at the m boundary (m £n). To evaluate 
F*(0) along path s', that starts on the side of the cylinder (Fig. 11*7), 
the directed radiative- energy intensity is needed and is written as 





m+i-1 



(11-38) 


(11-39) 


If this expression is integrated over the side of the cylinder, one 
obtains 

2tt b//b^+d2 

F^gCo) = 4 y* dcf.y' d(cos0) I^g(e,(})) , (11-40) 

0 *"^0 
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where d = a - r and a is the radius of the plasma* When the path s' 
m ^ 


orginates on the top surface, different regions are denoted by different 
s' (Fig. II-8). 


cos 6 


*/T — T ' 

+ ST 

" *1 m 


(11-41) 


3,m 


where 




T^cos^ - sin^(j) 


m 


sin6 


(11-42) 


here ^^e directed radiative-energy intensity is: 


1 L -? j=m 


i=m 


(11-43) 




where k=m+nfcnis the number of boundaries crossed to define 

0. , , and 

3+1, m* 

b/cos0 j * h 

0 j = m 

dH 


- 


r cosd) 
m ^ 


2 2 
T sin^(|) 
m ^ 




sin0 


j m,k. 


(11-44) 


If this expression is integrated over the top of the cylinder, we obtain 


IT/ 2 ®j+l 


FvtCO) = 4 


s . /«/ 


d(cos0) I^t(0,4)) , 


(11-45) 


cos 0. 


V 


u2^ 2 

b +a 






% 


I 


I : 




Figure H-8. Angles 0j^n,, 0j+i^m are defined by the and + i 

surfaces on the top of the cylinder with the origin at the 
n,th cylinder. 
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Equations (11-40) and (11-45) add up to give Fv(0). 

To obtain the backward contribution to F^(0), it is simpler to 
write the angular dependence in terms of 0 and A since these are the 
quantities that define the boundary regions of interest. Hence 
can be written as 


tt/ 2 ^j+l,m 

E / "‘ / 

* A v' •' n 




F‘ (0) = 4 

US'- 


d(cos0) l'g(e,A) 


( 11 - 46 ) 


^j,m=0 ^j,m 


where d = a + r^^ and 


m+n-2k-l 




( 11 - 47 ) 


i=m 


Here k = number of cylinders that are crossed twice, and we define 

lJSI 


r 


r cosA 
m 


2 . 2 . 
r_ sin A 
m 


sin0 


j - m+n - 2k 


s. 

J 


■ 


r cosA 
m 


V 


2 2.2. 

r. - r sin A 

J. 


sin0 


3 = m 

j I* m 
j < m+n 


2k , 


k.. 


and 


^i ■ + n - 2k ♦ 1 . 


( 11 - 48 ) 


Equations (11-46) through (11-48) give the contribution to F"(0) from 
the side. To obtain the contribution from the top, the cylinder will be 
divided into angular regions as seen in Fig. (II -8) and defined by 
Eq. (II-3). Hence, 
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If/ 2 ^j+1 


F^tCO) = 4 


COS0. = 


“i 

■ •! 


e. 

j 


'b^+a^ 


I‘^(6,A) d(cos0) . 


(11-49) 


and 


■ I I 


i *» f 

If 

ii 




s' 


(11-50) 


where k = m + number of boundaries along path at angle 0^ . Also, in this 
case, 

b/cos0 j = k 






r cosA 
m 




2 . 2 . 
r sin A 
m 


sin0 


r cosA - 
m 


_1j. 


2 . 2 . 
r sin A 
m 


sin0 


and 


V. 


■Yj = Vn • 


j > n, j k 


j - n, j ?< m 


j = m 


(11-51) 


(11-52) 


The sum of Eqs. (11-46) and (11-49) give F^(0). Therefore the sum of Eqs. 
(11-40), (11-45), (11-46) and (11-49) give F^(0). 

The only place it is of interest to calculate the current is on the 
surface of the plasma. Here J^(0,n) = J^(0,fi) so the previous analysis 
for calculating the total radiative- energy flux in the backward region 


\£ 
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is correct when Eqs. (11-46) and (11-49) are multiplied by . When 
the current is calculated along the x axis 

n • = sin6 cos(J) , (11-53) 

A numerical integration is carried out to obtain the current in this manner. 

When the core emits black-body radiation but the outer region remains 
semitransparent, the preceding analysis is correct for the shaded area in 
Fig. (I I- 9) because the path- lengths through this region do not go through 
the black body emitting core. However, for path lengths that pass through 
the core region, the directed radiative-energy intensity in Eq. (11-47) 
and Eq. (11-50) will become 


B 2 Y(s,,i-s^) 


1(0, A) = .11 e ^ 

V j=m 


m 




i=2 


i Y U. ,-s.) , 

e -1 ri g'r 3+1 r » (11-54) 

j=m 


where 


I® = 2h ' 




(II-55) 


and 



r cosA 
m 


V 2 2 . 2 , 

r. - r sin A 
j m 

sin6 


(11-56) 


With the above form of the directed radiative-energy intensity for paths 
that transverse the black-body core, the previous analysis can be used to 
calculate the total radiative-energy intensity. 

Now that the radiative energy-flux has been calculated it can be 
substituted into the rate equations to re-evaluate the state densities 
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within the plasma. The radiative flux can then be recalculated and this 
process is repeated out until convergence is obtained. Basically the non- 
linear coupling of the total radiative-energy flux and the state densities 

makes this iteration procedure necessary. 

In this chapter the rate equation formalism was first developed. 

Then the expression for the directed radiative-energy intensity was 
derived. Finally, expressions for the total radiative-energy intensity in 
a cylindrical geometry were obtained. 
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CHAPTER III. CODE DESCRIPTION AND TESTING 

A. Discussion of Code 

1. General Description 

A computer code was developed to solve the complex, non-linear prob- 
lem described in the previous chapter. In this code the rate equations 
Eq. (II-ll) and the total radiative- energy intensity Eq. (11-31) are 
treated as separate problems. The rate equations are first solved to 
give the source function and the attenuation for the radiative-energy 
intensities, then the radiative intensities are substituted back into 
the rate equations to recalculate the state densities. This process is 
repeated until convergence for the radiative energy intensities is obtained. 
A flow chart for this procedure is given in Fig. (III-l). The logic of 
the code is divided into four major sections and each section will be 
individually discussed. 

The first section initializes the data to be used in solving the 
problem. The next two sections evaluate the state densities and radiative- 
energy intensities, respectively. Finally the solutions are checked for 
convergence to determine if the problem is solved. 

2. Sectional Discussion of the Code 

In initializing the problem in the first section, the plasma must be 
described by its constituents and geometry. Also the number of the itera- 
tion must be given to allow for radiative energy-fluxes to be used from 
previous calculations if convergence was not obtained. 

Parameters describing the atomic composition of the plasma and the 
geometry of the system are input in the first section of the computer code. 
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Figure III-l. Flow chart of principal sections of computer code. The 
excitated state densities and total radiative- energy 
intensities are calculated in DENGAL and RADCAL, respectively. 
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The atomic data consists of the terms to be substituted into the rate 
equations, e.g. Eq. (II-ll). Here the data can be global or local. 

Global data consists of factors that are uniform throughout the volume 
of plasma. These terms include the Einstein A and B coefficients, energy 
levels of the plasma components and the frequencies of the various trans- 
tions. Local data describe conditions that exist within the specific 
regions of the plasma. Here the collisional reaction rates and line 
widths are found. 

Next the geometry of the system must be described. Cylindrical 
geometry is used in the present analysis. Thus the half height of the 
cylinder as well as the number of regions and their radii must be input 
into the computer code. The final step in the initialization process 
is to read in the number of the iteration and, if necessary, the electron 
state densities from the previous iteration. 

The rate equations Eq. (II-ll) can be solved once the plasma pro- 
perties are entered into the program. This is described in section two. 

A flow chart of the subroutine involved, RATEl, is shown in Fig. (III-2). 
First the collisional-radiative terms are combined into the form of Eq. 
(11-15) to be used in the matrix formulation of the rate equations. Next 
the coefficients are scaled down to the order or unity and the system of 
equations is closed using Eq. (11-12). 

With the state densities evaluated in GAUSSl, the radiative -energy 
intensities can be calculated in RADCAL. Figure (III-3) shows the sub- 
routines contained in CGYPAC where the calculations in RADCAL are actually 
done. Each of the subroutines in Fig. (II1-3) computes the different 
representatives of the directed radiative -energy intensity within the 
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RATE! 


Figure III-2. 



Flow chart of subroutine RATEl. This is the subroutine 
in DENCAL in which the excited state densities within the 
plasma are calculated. 








CGYPAC 
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Figure III-3. Flow chart of subroutine CGYPAC. This is the subroutine 
in RADCAL in which the total radiative-energy intensities 

are calculated. 
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plasma and DFPS evaluates the directed intensity along paths that originate 
on the side of the cylinder and are in the forward area. The form of the 
directed intensity is that of Eq. (11-38). DFPT calculates the directed 
intensity from the top of the cylinder in the forward volume and the form 
of this intensity is that of Eq. (11-43). In the backward volume DFMS 
computes the directed flux from the side of the cylinder and its form is 
that of Eq. (11-47). Finally DFMT evaluates the directed intensity from 
the top of the cylinder in the backwards direction and was the form of 
Eq. (11-50). Each of these must be integrated over the quadrants where 
they apply. Thus an integrati.ng routine to do double integrals containing 
a functional dependence for the inner integration had to be developed. 

A ten-point Gaussian- quadrature routine is used to do both the inner and 
outer integrals. To ensure intermediate data is not lost between iterations, 
the results of CYGPAC and RATEl are written on disks. 

The radiative-energy intensities are chosen to determine if the solu- 
tions converge. Because of their complex dependence on the state 
densities, the radiative-energy intensities converge more slowly than the 
state densities. The criteria used in these calculations are that no more 
than one percent of the solutions can diverge by more than 15 percent and 
only ten percent of the solutions can diverge between five and fifteen per- 
cent. Once these criteria are met the calculation stops and final solutions 
are printed out. 

B. Tests of the Code 

Within the computer code described in the previous section two major 
calculations are performed. The first is the calculation of the state 
densities in RATEl and then the computation of the radiative-energy 
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intensities in CYGPAC. Thus each calculation section was internally 
checked to determine if they were functioning properly. Then the entire 
code which contains both of these subroutines was tested to decide if it 
functioned properly. 

RATEl was first successfully tested on a system of three eqiiations. 

Then to see if a system with multiple electron states could be treated, 

a low-pressure helium plasma was treated. ^ ^ Helium was chosen in order 

to take advantage of the availability of electron cross-sections and 

state lifetimes. At this low pressure, 10 torr, the radiative-energy 

intensities in the plasma do not contribute significantly to the rate 

equations. These calculations illustrated the importance of the metastable 

f291 

states as was observed in experiments.'- 

To test the internal consistancy of CYGPAC a transition was chosen with a 
frequency of 4 x 10^^ sec with an inverse lifetime of 4 x 10 sec , 
an upper-state density of 1.5 x 10^^/cm^, a line width of 1.6 x 10 
sec"^, and an attenuation of 0.3 cm“^ . The cylinder had a 
height and diameter of 20 cm in these calculations. The first check 
involved varying the number of regions within the cylinder to determine 
how this affected the calculated radiative-energy flux at the edge of 
the cylinder on the centerline. The results of these calculations are 
seen in Table (III-l) for 1 to 5 regions within the plasma. The solution 
changes only 1.6 percent with the variation in the number of regions. 

This difference arises because the mesh of the integrations over the 
volume of the cylinder becomes finer as the number of regions is increased, 
thus the solution becomes more accurate. If a large multi-region problem 
were to be treated using this code then a better integration routine over 


Number 

of 

Regions 

F^(0) evaluated 
at the surface 

1 

38.5 ergs/cm^ 

2 

/ 2 

38.5 ergs/cm 

3 

38.3 ergs/cm^ 

4 

2 

38.2 ergs/ cm 

5 

, 2 

38.0 ergs/cm 


Table III-l. Comparison of the 
Radiative-energy Flux at the 
Surface of a Cylinder For a 
Variable Number of Interior 
Regions. 
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the azimuthal angle would be necessary to obtain accuracy. For example, 
instead of a 10-point Gaussian- quadrature integration, a 96-point quadra- 
ture could be used. 

As a further check, these results were compared to results for in- 
scribed and circumscribed spheres. The latter should provide lower and 

upper bounds on the total radiative-energy intensity and represent cases 

/ 

where an analytic solution is possible. As shown in Appendix A, the 
total radiative-energy intensity for a sphere is: 

(III-l) 

where 

ANhvg(v) (III-2) 

* ■ 47T * 

and a is the radius of the sphere. 

The radiative-energy intensity of the surface of an inscribed sphere 
containing the same plasma as in Table (III-l) is 21.9 ergs/cm^ while a 
circumscribed sphere has a radiative-energy intensity of 178.3 ergs/cm . 
Thus, the result from the code is indeed bracketed by these values. The 
reason for the large difference in the results for the two spheres is 
that the attenuation along the path length changes for each of these 
cases. In the case of the inscribed sphere the maximum path length is 
20 cm, while the circumscribed sphere’s maximum path length is 28.3 cm 
compared to a maximum path length of 22.4 cm for the actual cylinder. The 
ratio of the radiative- energy intensities is approximately 1:2:18 for the 
inscribed sphere to the cylinder to the circumscribed sphere. If just 
the exponential attenuation, e^^, is considered for the maximum path 
lengths in the plasma this ratio is approximately 1:2:22. Thus, CYGPAC 
is functioning properly. 


F (0) = 
Y 


2'rrf 

Y 


2ya 


- 1 ) -1 
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The previous tests prove that the results from each of the major 
calculational subroutines are correct. Now the code with both of these 
subroutines must be examined. The model used in these tests is a two 
species plasma in which species one is a fojr-lievel atom and species two 
is a three- level atom with 3 and 2 optical transitions, respectively. 

This simple model allows the state densities of the first iteration to 
be pre-calculated by hand as an additional check for the code. In the 
initial test in this series a case where there were high absorptions for 
all the transitions in the plasma was investigated. The high absorptions 
guarantee the total radiative- energy intensities will be small. Hence, 
the solutions will rapidly converge. Table (III-2) shows the state 
densities and radiative-energy intensities for this case. A converged 
solution is obtained in two iterations and the total radiative-energy 
intensities are small. 

Table (III-3) presents the state densities and calculated radiative- 
energy intensities for a case where a high gain is assumed. In this case 
the solutions rapidly converge because the saturation radiation intensity 
Ig is realized that is: 





(III-3) 


where c is the speed of light, Cy and are the inverse collisional 
life-times of the upper and lower state, respectively; A is the transition 
probability and B is the Einstein B coefficient for the transition. 

The final test in this series is a case where there are non- saturated 
gains present in the plasma. Table (111-4) shows the state densities and 


State Densities 
#/cm^ 


Radiation Intensity 
ergs/cm^ 


Iteration 

1 


Iteration 

1 

2 


1 

2 

It 

1.37x10^ 

17 

1.37x10 


1.0x10^" 

17 

1.0x10 

17 

1.13x10 

1.13x10^^ 


1.0x10^^ 

1.0x10^^ 

1.0x10^^ 

1.0x10^^ 





8.20x10"^ 

8.20x10"'^ 


4.44x10"^ 

4.44x10'^ 

s.esxio'"^ 

8.68x10"'^ 


1.81x10"^ 

1.81x10"^ 

3.52x10"'^ 

3.52x10“'^ 



.. . 


Table III-2. State densities and radiative-energy intensities for a 
case with large absorption. 
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Iteration 


Iteration 


State 

Densities 

#/cm^ 


5.48x10 


13 


4.52x10 


1.60x10 


1.00x10 


TI 


14 


14 


Radiation 

Intensity 

ergs/cm*’ 


7.36x10 


13 


6.06x10 


T3 


1.27x10 


14 


9.86x10 


13 


1.97x10' 


2.66x10 


3.81x10 


1.83x10 


2.66x10 


2.15x10 


1 


7.28x10 


6.00x10 


1.26x10 


1.01x10 


13 


TI 


14 


14 


3.64x10 


1 


2.66x10 


3.56x10 


7.28x10 


6.00x10 


1.26x10 


13 


1.01x10 


13 


14 


14 


3.61x10 


2.66x10 


3.54x10 


4.88x10 


13 


4.99x10 


13 


4.99x10 


13 


6.11x10 


1.12x10 


13 


6.01x10 


13 


6.01x10 


13 


12 


1.10x10 


12 


l.lOiclO 


1.03x10 


8.39x10 


9.92x10 


8.39x10 


9.92x10 


8.39x10 


4.99x10 


13 


6.01x10 


13 


12 


1.10x10 


12 


9.92x10 


8.39x10 


Table III-4. State densities and radiative-energy intensities for atoms with non-saturated gains. 
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total radiative-energy fluxes for this case. Four iterations are required 
before an acceptable solution is obtained because now there are large 
radiative -energy intensities that are not saturated., 

C. Conclusion 

In this chapter the computer code developed from the analysis in 
Chapter II was discussed. Then the tests of the principal computational 
subroutines were presented showing that the calculations produce results 
that are correct. Finally, the tests of the entire code were discussed 
proving that the code would treat all possible attenuations that are found 


in a uranium plasma. 
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CHAPTER IV. URANIUM PLASMA MODEL 

A. Introduction 

The first step in formulating a plasma model is to identify the 
principle constituents of the plasma. For the plasma under consideration, 
this entails identifying the concentrations of the various ions of uranium 
present. Calculations have been done at the University of Illinois by 

( 27 ) 

C. Bathke to determine the fractional densities of various species in 
the plasma and the resulting electron flux. He calculated the constituents 
of the plasma using the Saha equation as applied by Krascella^'^^ to 
uranium. Results quoted in the present work for the electron flux and 
plasma constituents are taken from Bathke 's work at a pressure of one 
atmosphere. 

The atomic constituents of the plasma change when the plasma tempera- 
ture changes as is shown in Table (IV-1). The ionization fraction ranges 
from 10 at 5000°K to 0.5 at 8000° K. Doubly ionized uranium can be 
neglected as a source of radiation in these calculations because of its 
low density. Although the neutron flux dramatically alters the high- 
energy electron flux, it does not significantly change the densities of 
the plasma constituents. This is because at these temperatures thermal 
effects are the main cause of ionization rather than the high energy 
electrons. 

In the modeling of the plasma both global and local properties must 
be considered. Global properties are those properties that are constant 
for a given species throughout the plasma. Examples of these are the 
spontaneous decay lifetimes of the atomic states. Local properties are 
constant for a species only within a given region of the plasma. An 
example is the electron flux which changes in each region. 
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i 


T Temp N++ 

(1000®K) 


5 

1.39 X 10^® 

3.56 X 10^^ 

3.51 X 10^° 

— — 

6 

9.86 X 10^' 

17 

1.15 X 10 

2.78 X 10 

7 

5.77 X 10^^ 

2.32 X 10^^ 

6.56 X 10^^ 

8 

2.58 X 10^^ 

17 

3.26 X 10*' 

14 

7.18 X 10*^ 


Table IV-1. The atomic densities (#/cm^) of the various species 

of uranium present in the plasma for various temperatures 
at a pressure of 1 Atm. (27j n© is the density of 
neutral uranium. N+ and N++ are the densities of 
singly and doubly ionized uranium. 
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B. Radiation Data 

Electron energy levels must be identified in each of the radiating 
constituents to develop a model. This can be accomplished if the tran- 
sitions in both neutral, U®, and singly ionized, U"*', are identified. Then 
by describing the transitions by their inverse lifetimes (Einstein A 
coefficients), absorption and stimulated emission can be calculated. 

The Einstein A coefficient is sufficient to describe radiative 
processes within the plasma because absorption and stimulated emission 
can be described by the Einstein B coefficients. The ratio of the 
Einstein A coefficient to the Einstein B coefficient is 


A. • o 1. 3 

ij _ STrhv 

B. . ■ 3 

c 


(IV-1) 


where A., is the inverse lifetime, B. . is the Einstein B coefficient 
describint stimulated emission, h is Planck's constant, v is the frequency 
of the transition and c is the speed of light; also 


B. . g. 

= h. 

B. . g. 

ji *1 


(IV-2) 


where B^^ describes absorption and is the ratio of the lower to 

upper state degeneracies. 

The radiative term in the rate equation corresponding to absorption 
or stimulated emission is of the form 

PvCO) 


^ ^ ®ij c 


(IV-3) 


Thus once the inverse lifetimes are known, the radiative processes within 
the plasma can be described. 
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The first reported data^^^^ describing transitions in uranium was 
from an experiment at the National Bureau of Standards (N.B.S.) where the 
plasma temperature was 5000®K. Intensities for two hundred and eighty- 
five transitions in U® and two hundred and seventy-eight transitions in 
U"*" were reported. These workers did not, however, report lifetimes. 

Thus it is necessary to estimate A coefficients as follows. From the 
measured intensities, relative inverse lifetimes (i.e.,A^/A 2 ) can be 
obtained using the relation: 


I2 ’ 


(IV-4) 


where I is the relative intensity, v is the transition frequency, N 
is the atomic-state density of the emitter and 1 and 2 refer to transition 
1 and 2. Since the atoms were in thermodynamic equilibrium, Boltzxiann 
statistics is assumed valid, in which case N 2 /N^^ becomes: 


^2 -E/kT 

NT" ® 


(IV-5) 


where E is the difference in energy between state 2 and 1. 

In uranium the transitions are governed by j-j selection rules as 
opposed to 1-s selection rules in hydrogen. Figure (IV-1) shows the 
range^^^^of known levels in U® and U^. In neutral uraniinn a single family 
of upper states is identified that has even parity. These states emit 
radiation to two families with lower energies tliAt have odd parity. The 
same situation exists in singly ionized uranium but here there are three 
lower energy families with odd parity and one high-energy family with odd 
parity that can emit radiation to the lower lying even parity states 
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4 4 - 1 

(e.g. £ p **• f d). In this figure the state energies were given in cm , 

the conversion between cm”^ and eV being 

1 eV = 8066.1 cin"^ (IV-6) 

Later experiments at Los Alamos^^^'^^^ provided additional data on 
other transitions in uranium. Although many more transitions were identi- 
fied in these works, most of them could not be used because data on the 
conditions of the plasma were not described, (i.e.^ equilibrium vs non- 
equilibrium) . Consequently, only transitions that originated from the 
same upper state, making N 2 /Nj^ = 1 in Eq. (IV-4) could be added to the 
existing A values with any confidence. This process added 982 transitions 
to U® and 298 transitions to U^. The data are summarized in Appendix B 
(Tables B-1 and B-4) . 

(32) 

Thus the N.B.S. experiment identified a set of basic transitions 
in uranium and the Los Alamos experiments provided data for 
additional transitions. An example of this identification process can 
be seen in Table (IV-2) for some transitions in U*. Here the number of 
transitions positively identified by their lifetimes increased by a value 
of 370% using the LASL data, but some A values are still uncertain. This 
inability to include all the possible transitions is one of the major 
shortcomings of the present work. Still, the omission does not seem to 
be overly crucial, since the intensities of the missing lines were not 
predominant in the experimental observations. 

Absolute A values must be known for the calculations. To obtain 
these from the ratios, one of the lifetimes in each atom must be known. 

For U° the relative inverse lifetimes are converted to an absolute scale 

- 1 8-1 

using the 27887 to 0 cm absolute A value of 1.37x10 sec obtained 






FINAL STATE 


I 28444 

N 

I 28454 

T 

28470 
28504 
28543 
28563 
28650 

Table IV-2. Relative inverse lifetimes for some states of UI based on Ref. 32, Ref. 33 and 
Ref. 34. Energies are in cm~^. 
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by Klose,^^^^ e.g., this gives a value of 1.30x10^ sec ^ for the 28504 
to 620 cm”^ transition. No absolute A value was available to normalize 
the relative A values in U'*'. There is a state coincidence at 27499 cm ^ 
in U° and U^; so assuming a 50:50 ratio of U°:ll* Nj^/N2 will become 1 
in Eq. (IV-4) if both transitions originate from these states. This 
evaluation of A values from the 27499 cm ^ state in U establishes an 
absolute scale for A values in ionized urani\mi. 

In the above manner the state energy levels were identified and the 
Einstein A coefficient for the various transitions was calculated. A 
summary of the lifetimes used is presented in Appendix B. 

C. Collisional Data 

In addition to the radiative data, the collisional cross sections 
are needed to describe the atomic states in the plasma. Atom-atom and 
atom-ion collisions can be neglected because the atoms are at thermal 
energies and due to their large mass are moving slowly compared to the 
orbital speed of bound electrons. Hence the electrons will adjust them- 
selves to the changing force field in collisions with atoms without 
undergoing transitions.'- ^ 

1. Ionization Cross-Sections 

Direct experimental measurements of cross-sections for electron 
excitation and ionization in uranium are not available; thus, the necessary 
data had to be calculated. The ionization and excitation cross-sections 
by electron impact have been calculated from formulae proposed by 
Vriens^^^^ based upon the Gryzinski model collected for enhancement due 
to ions^^^^. Vriens corrects for the asymmetry in the original Gryzinski 
cross-sections by allowing the incident electron to gain kinetic energy 
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equal to the potential energy it loses as it interacts with the atomic 

(36) . 

electron. The resulting formula for the ionization cross section is 


■n' 

■ E + 2UI 


‘A A . 2 A Ui\ ME/m.1 

\UI " eJ 3 \UI g2/ ■ E + UI J 


(IV-7) 


where UI is the ionization potential of the state under consideration. In 

the present work, however, the factor (E + 2UI) ^ in Eq. (IV-7) is replaced 

by e"^ as well as in the excitation cross-section formula presented later 

(37) 

in order to account for the focusing effect of ions. 

The accuracy of these cross-sections compared to experimental cross- 

sections is found to be about a factor of two. The excitation cross- 

sections presented in the next section are found to be accurate to an 

order of magnitude. The large discrepancy in this case is due to the 

smaller energy differences possible in excitation events as opposed to 

ionization events. 

Figure (IV-2) shows the calculated ionization cross-sections for the 
ground state of neutral and singly ionized uranium. The cross-section for 
neutral uranium peaks earlier and is larger than the corresponding cross- 
section for singly ionized uranium because of the larger ionization potential 
of singly ionized uranium. These characteristics are again observed for 
the cross-sections found for the higher energy states in neutral uranium 
shown in Fig. (IV- 3). 

2. Excitation Cross-sections 

The excitation cross-section is more complicated in the Gryzinski 
model, since the excitation to the final state and the next higher energy 
level in addition to the ionization potential are necessary to describe 


CROSS-SECTION, cm 
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the interaction. The resulting cross-section can be written as 



where is the excitation energy to state n, is the excitation energy 

to state n+1 and R = 13.58eV. The rise from threshold to peak value is 
determined by the energy of the next higher state and over this energy 
range, £ E £ is replaced by E in Eq. (IV-8). 

Because of the smaller energy transfers involved in excitation 
collisions, these cross-sections will generally be larger than for ion- 
ization. Also, as can be seen from Fig. (IV-2) the rise in the ionization 
cross-section from threshold to peak value takes place over an energy 
range of the order of electron volts. In contrast the rise from threshold 
to peak energy for the excitation cross-sections takes place over a smaller 
energy corresponding to the difference in energy between the n^^ and the 
n^^+1 state. This sharp rise is evident from Fig. (IV-4) and Fig. (IV-5). 

Figure (IV-4) shows cross-sections for excitation from the ground 
state to three different higher energy states in neutral uranium. The 
largest cross-section has the lowest excitation potential and the peak 
values of the cross-section decreases as the excitation potential increases. 
These trends are seen in Fig. (IV-5) which shows the cross-sections for 
various electron excitations to the 3.9959 eV energy level in neutral 
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Figure IV-4 . Various excitation cross-sections from the ground state 
of neutral uranium. 
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Figure IV-5 . Excitation cross-sections for various states to the 
3.9192eV state in neutral uranium. 
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uranium. The increase in peak value of the cross-section due to the 
decrease in excitation potential is further enhanced by the decrease in 
ionization potential as seen by comparing the excitation cross-section 
from the ground state to the .5301 eV state in Fig. (IV-43 to the exci- 
tation cross-section from the 1.6792 eV state in Fig. (IV-S), Both 
of these cross-sections have approximately the same peak value although 
the excitation potential is doubled in the latter case. The difference 
in the ionization potential of the initial state compensates for the 
difference in excitation potential, with the atomic parameters now 
known, a model for the plasma can be implemented. 

D. State Reductions 

Due to limitations on the size of the computer program, all states 
in uranium cannot be treated explicitly. Hence^ the number of states to 
be considered were reduced to a smaller number of "composite" states. 

To understand this reduction, consider the rate equation for state p 

E (r +c 3 n + V CR +c 3 n = (iv-93 

pq pq P qp qp q P » 

n q 

where q represents other states in the plasma, R the radiative terms in 
the rate equation, C the collisional terms in the rate equation, D the 
recombinational terms from the rate equation and N the atomic state 
densities. 

For the model of uranium to accurately represent the actual plasma, 
the coefficients of the composite states in the rate equations must be 
equivalent to the coefficients of their constituent states. In the 
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present model, the composite transitions represent the average of all 
transitions from the constituent states. Hence, the A value and cross- 
sections of the composite state must be comparable to those of the 
constitutent states. In the present work, the composite state lifetime 
did not vary by more than a factor of two from the lifetimes of the 
constituent states. Also, the energy of the constituent states was within 
200 cm ^ of the energy of the composite state. This ensures the cross- 
sections of the composite state will be equivalent to the cross-sections 
of the constitutent states. Figure (IV-5) illustrates this because the 
smallest two cross-sections are for states that are separated by 619 cm 
Figure (IV-6) shows a two level diagram where both the upper and lower 
levels represent composite states. 

Using the above criteria the number of states was reduced from 169 
to 62 in neutral uranium and from 95 to 46 in singly ionized uranium 
while still treating 822 optical transitions explicitly. 

E . Elect r on Flux 

In the preceding sections of this chapter global physical parameters 

were discussed. In this section the electron flux which determines local 

conditions in this analysis will be discussed. The calculations discussed 

f27') 

in this section are taken from the work of C. Bathke 
1. Temperature Dependence 

The first parameter we study is the temperature dependence of the 

electron flux for a given neutron flux. The reference neutron flux level 

14 2 

used in this analysis is 2x10 neutrons/cm -sec. The range of tempera- 
tures of interest is between 5000®K and 8000®K and the pressure is one 
atmosphere. 


1 



composite model actual states 

Figure IV-6. Schematic diagram of two level composite state 
representing five constituent states. 



Figure (IV-7) shows the variation in electron flux with different 
temperatures. The non-Maxwe Ilian component of the electron flux decreases 
as the temperature increases due to the corresponding decrease in the 
fission cross-section over this temperature range. The reduced fission 
cross-section causes the source of non-Maxwe Ilian electrons to decrease, 
resulting in approximately an order of magnitude reduced in the high- 
energy electron-flux in going from 5000*K to 8000°K. This, coupled with 
the lower densities in the Maxwellian flux at low temperatures, implies 
that more radical non-Maxwellian excitation will occur in the low- temp- 
erature case. 

2. Neutron Flux Dependence 

A much mo3,‘e dramatic variation in the electron flux is seen when the 
exciting neutron flux is varied. In this study the range of variation is 
from 10^^ to 10^^ neutrons/cm^-sec. Since the reaction rate for fission 
is proportional to the number of neutrons in the system, the change in 
the electron flux closely follows that of the neutron flux as can be 
seen in Fig. (IV-8) . 

The electron flux is so important because within the plasma all 
collisional reaction rates are calculated using it. The parametric study 
done here is of interest because it will be a direct test of the 
effectiveness of the non-Maxwellian tail of the electron flux in 
exciting the atoms in the plasma. 

F . Data Processing 

The preceding sections of this chapter have dealt with the various 
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Figure IV-7. Electron flux for a neutron flux of 10 n/cnr-sec at 
temperatures of 5000°K, 6000°K, 7000“K and 80n0°K.^^^^ 
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physical phenomena that describe the plasma. In the present section the 
synthesis of these processes will be discussed. 

Two types of coefficients are seen in the rate equation [Eq. II-ll)]; 
those that vary with and those that remain constant with respect to the 
radiative- energy intensity. Collisional terms in the rate equations do not 
change with the radiative-energy intensity because the variable here is 
the electron flux. Thus within a given region these processes need only 
be calculated once. Throughout the plasma the state lifetimes for 
spontaneous decay, Einstein A coefficients also remain constant. Hence, 
these terms can be combined and saved for any calculations within a given 
region. This saves time when iterations are done in the code. 

Absorption and stimulated emission depend upon the radiative-energy 
intensity and have to be calculated each iteration. These terms, combined 
with the constant terms, provide all of the coefficients employed in the 
rate equations . 

On:e the state densities are evaluated from the rate equations, the 
only parameter required to calculate the radiative-energy intensity is the 
line width. Due to the high pressure of the plasma a combined Doppler 
and collisionally broadened line widths are used^‘^^’'^^\ Thus the line 
width, FWHM, is written as: 




Ao = 


+ Av , , 
V coll 


ClV-iO) 


where V is the transition frequency, is the collisionally broadened 

width and AVj-j is the Doppler width given by 


AVp = 


V 


8kT In 2 
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where k is Boltzmann's constant, T is the temperature of the plasma and 
M is the mass of the atom. The collisional width used in this work was 
developed by Linholm^"^^’^^^ , who, by including the time dependence of the 
collision into the interruption broadening theory of collisions, accounted 
for close collisions. His results reduce to the interruption broadening 
results near line center or at low pressures. Also at high pressures or 
at frequencies in the wings, his result reduces to the statistical 
broadening result appropriate for these limits. Thus his model correctly 
describe the line shape in these limits. The Linholm collisional width is 


^^coll = ^ N » 
o 


(IV-12) 


where e is the electron charge, m is the electron mass, f is the oscillator 

strength, and N is the density of the atomic species. The oscillator 
strength can be written 


Ac^ m 
STT^v^ 


CIV-13) 


where A is the transition probability. By substituting Eq. (IV- 13) into 
(IV-12) , the collisional width becomes 


Ac^N 


(IV- 14) 


To obtain the total line width Eq. (IV-11) and (IV-14) are substituted 
into (IV- 10) to give 


Av = -§ILMLi£2 V + Ac^N 
Mc^ 16 ttv^ 


(IV-15) 
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The ratio of the collisional to Doppler contributions to the line width 
is obtained by dividing the first term in Eq, CIV- 15) into the second 

^^coll _ / M Ac^N ^ (IV-16) 

Ayp ^ 2 it1cT 32ttv^ 

As can be seen here when the frequency of the transition is small or the 
density or transition probability is high, the collisional width becomes 
more predominant. Conversely, when the density or transition probability 
is low or the frequency is high then the Doppler width predominates this 
expression. Equation (IV-15) is the final parameter necessary to cal- 
culate the radiative- energy intensity. 

G. Conclusion 

In this chapter tlie distinction was made between global and local 
properties. Then the global properties were explained and used to give 
a model for the uranium plasma. This model is based on the modified 
Gryzinski cross-sections and experimentally measured radiation intensities 
Thus, the model is limited by the accuracy of the theoretical cross- 
sections and the calculated lifetimes. In the absence of other data to 
incorporate into a model of uranium, this model is the most complete 
to date. Next, the electron flux, which defines the local plasma 
properties, was discussed. Finally, the combination of these properties 
in the transport code was explained. 
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CHAPTER V. RESULTS AND DISCUSSION 

A. Physical Parameters and Convergence 

Chapter II described the formalism that was incorporated into the 
computer code discussed in Chapter III. This chapter presents the results 
of these calculations. Unless otherwise specified, the plasma is at a 
xiniform temperature and has a pressure of one atmosphere. The uranium 
plasma, consistent with projected gaseous core experiments, is cylindrical 
with a diameter and height both equal to one meter. 

Before these results are presented, the convergence criteria used in 
the code must be discussed to understand the validity of the calculations. 

A choice of either the atomic state densities or the radiative-energy 
intensities could be used to determine convergence. If there is no 
radiation incident upon the surface of the cylinder, the directed radiative- 
energy intensity can be written 
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I(v,oj = exp 

1 y(V,S* jds' 

1 f(v,s')ds* exp 
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(V-lj 

The source term, f(v,s) is proportional to the upper state density 
and the attenuation, Y(v,s) is proportional to the difference in the 
upper and lower state densities. Due to this complex dependence, the 
radiative-energy intensities are quite sensitive to changes in state 
densities and thus are used to establish convergence. 
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If the plasma is homogeneous so the source term and the attenuation 
are independent of position, Eq. (V-l) becomes 


ICv,o) 


Y(v) 


- 1 ). 


(V-2) 


The largest path length in the plasma is 141 cm. Assuming y changes 
from -.05 to +.05 between iterations, the exponential, evaluated for the 
maximum path length in the plasma, in Eq. CV-2) changes from .0009 to 1153. 
Due to this large change, two convergence criteria were used. The first 
criterion is then no more than one percent, i.e., 8 of the transitions 
can diverge by more than 15%. This accounts for transitions that are 
changing from small absorptions to small gains. The second criterion 
is that no more than three percent, i.e., 25 of the transitions, can 
diverge between 5% and 15%. Wlien both of these conditions are satisfied 
the solutions are assumed to be converged. 

Details of the convergence for the case of a flux of 2x10^^ 

2 

neutrons/ (cm -sec) and a uniform plasma temperature of 7000°K are 
shown in Table (V-l) . Column one gives the number of the iteration] 
the second column, N15, gives the number of transitions that have 
diverged between 5-g and 15%; and the third column, N15, gives the number 
of transitions that have diverged by more than 15%. The final column, 

N, represents the number of transitions that have converged to better 
than 5%. Initially, the convergence is slow but after the tenth iteration 
rate increases rapidly as the solutions stabilize. Finally, the rate 
decreases after a majority of the solutions have converged because only 
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a few transitions have yet to converge. By the eighteenth iteration 
96.5% of the transitions have converged to better than 5% in this case. 
Thus the observables calculated here should also be good to 5%. 


This method was used in every calculation and Table (V-2) shows 
the final convergence for the uranium plasma when the temperature 
varied from BOOO^K to 8000“K with the flux fixed at 2x10^^ neutrons/ 

(cm^-sec) . The 5000‘*K case converges so well because there are no 
inverted transitions, hence>the total radiative-energy intensities 
are small and do not make significant contributions to the rate 
equations. This effectively linearizes the problem because the 
atomic state densities are no longer dependent upon the radiation within 
the gas. As the temperature increases, the radiative-energy intensities 
increase, causing significant contributions to the rate equations and 
making the problem non-linear. Hence, the convergence becomes more 
difficult. In all cases the final solutions have converged to better 
than 5% for over 96.5% of the transitions in the plasma. 

B. Radiative-Energy Current 

In this section the radiative-energy current from a uranium plasma 
with a neutron flux of 2x10^^^ neutrons/ (cm^- sec) and a temperature variation 
of 5000“K to 8000®K will be discussed. Using the calculated atomic state 
densities, the radiative- energy current is computed from Eq. (11-29) where 

• fl' = sin0 cos4> , 


(V-3) 
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and the current is evaluated at the mid-plane of the cylinder. The 
current is evaluated at a frequency corresponding to line center of each 
transition. 

The Planck function. 



(V-4) 


describes the radiative-energy current from a system in thermal equi- 
librium. Equating Eq. (V-4) With the calculated radiative-energy 
current, the line "temperature", Tj^, can be written: 



(V-5) 


where J is the calculated radiative -energy current. T. when substituted 
into the Planck function will give the calculated radiative -energy 
current for the line at line center. Once the line temperature is 
calculated for all lines within the plasma, an effective temperature 
is defined by 


T 


eff 


<Tj^> . 


(V-6) 


Figure (V-1) illustrates this process. The solid curves represent the 
calculated radiative-energy current, J, from which the line temperature 
was calculated. These line temperatures are averaged to give the effective 
temperature. The current corresponding to the T^^^ defined as J 

eff 
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represents a weighted average as illustrated in the figure. T^££ is 
useful since the radiative -energy current from lines not included in 
this analysis can be computed by substituting the effective temperature 
and the frequency of the transition into Eq. (V-4) . This is illustrated 
by the dashed curve in Fig. (V-1). Finally to solid line is the 
Planck function evaluated at T 

eff 

Figures (V-2) through Fig. (V-5) show the J's plotted with respect 
to the Planck function at for a uranium plasma with a temperature 

variation of 5000°K to SOOG^K and a flux of 2x10^"^ neutrons/ (cm^- sec) . 
The horizontal curve is the Planck function at T^££ and the vertical 
lines are drawn from the Planck function to the radiative-energy current 
at line center as illustrated in Fig. (V-1). For 5000“K, no inversions 
are found within the plasma and T^^^ is 9108®K. Much hotter effective 
temperatures are found at higher plasma temperatures, as in the 6000" K 
case where T^^£ is 36,960"K. Along with the significant rise in 
a number of inversions are found in this case [see Fig. (Y~3)]. The 
most severe elevation in T^^^ occurs at a plasma temperature of 7000" K 
where T^^^ becomes 389,900"K [see Fig. (V-4)], Finally, a T^^^ of 

180,800"K is calculated for the 8000°K plasma case [see Fig. (V-5)]. 

14 2 

The T^££ variation for a flux of 2x10 neutrons/ (cm -sec) is 
shown in Fig. (V-6) where T^££/T is plotted against T, the plasma 
temperature. To show that this variation can be predicted and explained, 
the simple three-state model illustrated in Fig. (II-2) is considered. 
The 3.9605-eV to .0769-eV transition in neutral uranium is assumed. 

In a homogeneous plasma radiative-energy current is 
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Figure V-2. Radiative- energy current versus frequency 
measured with respect to the Planck current 
at an effective temperature of 9108®K. 
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Figure V-3. Radiative-energy current versus frequency 
measured with respect to the Planck current 
at an effective temperature of 36,960®K. 
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J (x) = ft • ^ (1 " ® 
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(V-7) 


If Y (y)x is large then the radiative-energy current will be proportional 
L 

to £ (vVy (v)* Substituting this expression into Eq. (V-4) the line 
L b 

temperature becomes 



where is the upper state density and is the lower state density. 

Table (V-35 suraiarites the state density and line temperature dependence 
on the plasma temperature for this case. This single transition will not 
exactly exhibit the variation found earlier for T^jj because the latter is 
an averaged quantity, but the line temperature shows similar trends. 

There is a rise in the line temperature as the plasma temperature varies 
from 5000“ K to 7000° K and then the line temperature decreases as the 
plasma temperature varies to 8000°K. If the system were in thermal 
equilibrium, the line temperature would be equal to the plasma temperature. 
This can be seen by substituting the Boltzmann factor, e ^ , for 

in Eq. (V-8). 

Two processes in the plasma cause the peaking effect observed in 

T The first is recombinational excitation within the plasma. In 

eff 

Table (V-3) the number of electrons in the upper state increases when 
the plasma temperature charges from 5000°K to 6000°K (since the non- 
Maxwellian tail in the electron distribution is inversely proportional 
to temperature [see Fig. (IV-7)], the non- equilibrium electron excitation 


t/k 

M -3 

Ny,cm 

Nj^tcm"^ 


Teff/T^ 

5000 

1.508x10^^ 

3.173x10^^ 

8415 

1.683 

6000 

4.793x10^^ 

1.843x10^^ 

12334 

2 . 056 

7000 

4.399x10^^ 

9.602x10^^ 

14598 

2.085 

8000 

9.36x10^"^ 

3.969x10^^ 

12011 

1.051 


TABLE V-3. Upper and lower state density and 

14 

line temperature for a flux of 2x10 
neutrons/ Ccm^-sec) calculated from a 
three state model. 




87 


is not caused by the non-Maxwellian tail, but by recomb inational 
excitation) . The second effect is the changing constituent densities 
within the plasma [see Table (VI-1)]. As the temperature increases, 
the number of ions and electrons in the plasma also increase causing 
more recombinational excitations to occur. If the number of neutrals 
remained constant or increased with temperature, the effective temp- 
erature would monotonically increase. The number of neutral atoms 
decreases as the plasma temperature increases, thus, the number of 
electrons in the excited state actually decreases [see Table (V-3)], 
as the temperature increases from 7000“K to 8000®K causing a drop 
in the effective temperature. 

The total radiative-energy current carried in the lines calculated 
over frequency intervals of 1x10^^ Hertz, is summarized in Table (V-4) . 

A typical lasing line width of 1x10 Hertz was used for lasing transitions 
in these calculations. Low frequency transitions are predominantly 
found in neutral uranium while high frequency transitions dominate in 
singly ionized uranium. 

The transitions in the 1 to 2 x 10^^ Hertz interval are all from 

14 

transitions in neutral uranium and the transitions in the 11 to 12 x 10 

Hertz interval are from singly ionized uranium. The density effects can 

14 

be seen directly in those two intervals. In the 1 to 2 x 10 Hertz 
interval, the energy carried in the transitions decreases as the 
temperature increases. This is because the neutral uranium density 
decreases as the temperature increases [see Table (VI-I)]. The energy 
carried by transitions in the 11 to 12 x 10 Hertz range increases with 


7 -8 


8 -9 


9 -10 


10 -11 


11 -12 


5000 

6000 

7000 

8000 

' 2.27x10^ 

2.38x10^ 

1.49x10® 

8.76x10® 

1.20x10® 

1.44x10® 

1.35x10^ 

7.58x10® 


2.31x10^ 

1.30x10^ 

6.78x10® 

2.14x10^ 

1.92x10^ 

2.42x10^ 

2.65x10^ 

5.02x10^ 

2.87x10^ 

5.28x10^ 

7.13x10^ 

1.62x10^ 

1.49x10^ 

3.18x10^ 

3.92x10^ 

2.15x10^ 


9.10x10^ 

1.40x10® 

1.07x10^ 

2.71x10^ 

5.49x10^ 

7.64x10^ 

4.28x10® 

2.20x10^ 

8.78x10® 

1.87x10^ 

6.56x10"^ 

1.00x10^ 

6.47x10® 

. 

9.06x10® 

1.19x10^ 

2.26x10® 

1.76x10® 

4.42x10® 


TABLE V-4. Radiative- energy current 
carried in the lines 
within the plasma with 
the units of ergs/cm . 
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temperature because the density of singly ionized uranium increases with 
temperature. Within a particular frequency interval, the total radiative 
energy current depends upon the number of linfcs in the interval and 
whether the transitions originate in neutral or singly ionized uranium. 


C. Opacity Calculations 

In this section the opacity of the plasma is calculated. This 
quantity is a measure of the absorption within the plasma and can be 
evaluated from 


V 



(V-9) 


where p is the plasma density and y is the absorption coefficient as 
given by Eq. (11-19) fot the lines and Eq. (11-21) for the continuum 
between the lines. The opacity is essentially the mass absorption 
coefficient of the plasma. Thus, at a given density, the larger the 
opacity the less radiation will penetrate the plasma. Likewise as the 
opacity increases, the emitted radiation will come from a smaller 
volume of plasma near the surface. 

All the opacities are calculated and then averaged over 1.55 x 
13 

10 Hertz frequency intervals according to: 


K Av 

<K > = - 4 — 
V Av 


(V-10) 


where Av is the line width and the sum is over all non-inverted lines 
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in the frequency interval plus the continuum between the lines. 

The dots in Fig. (V-7) are the result of such an averaging for a flux 
of 2x10^^ neutrons/ (cm^- sec) and a plasma temperature of 5000®K. 

Parks calculated the opacity for a uranium plasma at a tempera- 
ture of 5100®K and pressure of 3.58x10"^ atmospheres using experimental 

(32,44) 

intensity data from the National Bureau of Standards . The 

dashes in Fig. (V-7) represent his results corrected for the difference 
in densities of U“ and U^. 

In his calculations there was a 40/60 mixture of neutral to 

singly ionized uranivira while in the present calculations there is a 

97.5/2.5 mixture. At low frequencies the transitions are principally 

found in U° while at high frequencies the transitions are mainly found 

in U^. Since the opacity is the absorption divided by the density 

[Eq. (V-9)], Parks' opacities at high frequencies were divided by 24 

and at low frequencies they were multiplied by 2.44 to correct for the 

difference in the mixtures. The cut-off between high and low frequencies 

14 

was taken at the largest opacity at 8x10 Hertz. 

Figure ,(V-8) through Fig. (V-10) show the opacities for temperatures 
other than 5000“ K. The effect of the changing constituent densities can 
be seen in these figures. As the temperature increases the U“ density 
decreases. Consequently, the opacity decreases at low frequencies with 
increasing temperature. Similarly, the opacity at high frequencies 
increases as the temperature increases from 5000“ K to 7000“K due to an 
increase in the density of u"^. Also a dip in the opacity at approximately 







OPACITY, cmVgm 



Figure V-8. Average opacity ver^s frequency for a plasma at a 
temperature of 6000* K. 



OPACITY, cmVgm 



Figure V-9. Average opacity versus frequency for a plasma at a 
temperature of 7000*K. 


(O 


OPACITY , cm2/gm 



FREQUENCY, lo'^ hertz 




Figure V-10. Average opacity versus frequency for a plasma at a 
tenroerature of 8000®K. 
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l.xlO^^ Hertz that is observed in Parks calculations [see Fig. (V-7)] 
is also seen in the present 6000®K and 7000*K cases. The decrease in 
the opacity at high frequencies in the 8000®K case. Fig (V-10) is 
due to the increased excitation in singly ionized uranium which causes 
the magnitude of the absorption due to the lines to decrease [see Eq. 
(11-20)]. This is because the absorption is proportional to the 
differences between the lower and upper state densities. 

D. Flux and Size Considerations 

In the uranium plasma, the non- equilibrium effects are mainly due to 
recombination in the plasma. Hence, a change in the non-Maxwellian tail 
of the electron distribution function has no significant effect on the 
state densities in the plasma. Consequently, when the flux was varied 
from 2x10^^ neutrons/ (cm^-sec) to 2x10^^ neutrons/ (cm^- sec) at a 
plasma temperature of 6000®K, no significant changes in the radiative- 
energy current or the atomic state densities in the plasma occurred. 

Thus, the radiative-energy current calculated for one neutron flux in 
uranium is adequate for a range of neutron fluxes. Also, because the 
optical depth is small, less than a centimeter, the radiative-energy 
current does not change with a change in size of the plasma. Thus the 
total radiative-energy radiated by the plasma will be directly proportional 
to the surface area of the plasma. Then the radiative-energy current 
calculated here need only be multiplied by the surface area to obtain the 
total radiative-energy. 
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In this chapter the physical dimensions of the plasma were described 
and then the convergence criteria for the code were discussed. Next the 
results of the radiative-energy current calculations were presented, first 
introducing the concept of an effective temperature and then the actual 
energy carried in a frequency interval. Then the opacities of the plasma 
were presented along with a comparison to Parks ^ ^ earlier results for 
a 5000“K plasma. Finally the effects of changing the size of the 
plasma and the neutron flux were discussed. 


CHAPTER VI. CONCLUSIONS 


A. Uranitmi Calpulations 

The results of this work show that non-equilibrium effects in the 
plasma cause larger radiative- energy currents at the lines than would 
normally be expected. United Aircraft 's^^^ assumption of a black-body 
spectrum is true for the interior of the plasma but there is an important 
contribution to the spectra due to bound-bound transitions in a sheath 
near the surface. The energy carried in these transitions is about 2% 
of energy carried by the black-body spectrum in thxs frequency 
range. For calculating bulk quantities within the plasma a black-body 
spectra is sufficient; but if radiation processing or other frequency 
dependent processes are of interest, detailed analysis of the lines is 
necessary to correctly predict the emitted radiation. 

C44') 

In recent studies at United Technologies, Krascella'' found that 

there was up to an order to magnitude difference between the integrated 

(45) 

intensity and the black- body intensity at some temperatures, and Roman 

found the ground s_ate density in neutral uranium was depleted to a 

greater extent than was predicted by equilibrium calculations. He found 

the ground state density to be 8x10 atoms/cm and 2x10 atoms/cm in 

1/1 

two instances compared to 10 atoms/cm predicted by equilibrium 
calculations. This depletion of the ground state is indicative of non- 
equilibrium excitation within the plasma as illustrated by the non- 
equilibrium radiation levels observed in these calculations. 
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Rogers et predict a rise in the black-body radiative 

temperature corresponding to the outward heat flux above the gas 
temperature for the plasma core reactors they examined. Here a much 
greater rise in is observed than Rogers reported because he 
considered bluk quantities while the present work concentrates on 
line radiation where non- equilibrium effects would readily appear. 

At high temperatures, 6000*K and above, approximately 20% of the 
transitions were inverted. These inversions were essentially due to 
recomb inational effects. Due to the limitations of the model, exact 
inversions cannot be predicted. However, with the large ion and 
electron densities in the plasma, recombinationally driven inversions 
should be expected. Also observed was the effect of two constituents 
within the plasma. This is readily seen in the opacity calculations 
where the high frequency results were dominated by and the low 
frequency results were dominated by U®.; 

Finally due to the density of the plasma, line radiation comes from 
a narrow sheath on the plasma surface. This effect was assumed by 
Rogers and Latham in their modeling where they allowed a 1.4 cm 
sheath that would account for all optical radiation. In the present 
calculations, 1.4 cm would be greater than a mean free path for any 
of the line transitions in the plasma. 

The accuracy of these calculations depends primarily on the accuracy 
of the uranium model. Gryzinski^^^^ cross sections are used, but they 
have been successfully employed by other researchers previously (e.g.. 


used the model for helium). However, comparisons of the calcu- 
lated and experimental cross sections for cesium, which is similar to 
uranium in electronic structure, reveals that the calculated cross- 
section may be somewhat low. The uncertainty is estimated to be a 
factor of 2. This suggests that experimental measurements are badly 
needed. Also, the lifetimes of the atomic states should be measured 

directly to improve the accuracy of the results. Some progress is 

, . u. (33,50) 

being made in this direction 

B. Implications to Uranium Plasma Program 

The energy in line radiation represents a pertubation to the energy 
carried in the black-body spectrum of the plasma. Although only a 
pertubation, the radiative- energy current at a specific frequency corres 
ponding to a line center can be orders of magnitudes larger than the 
current expected from equilibrium calculations. 

In an uranium plasma, electron excitation may result from both 
direct excitation and recombination. Recombination will be the pre- 
dominant cause of inversions in a plasma where the constituent atoms 
have small excitation and ionization potentials. This is because the 
Maxwellian electron flux will cause excitation from ground due to the 
small energy differences and non- equilibrium excitation will be caused 
by recombining electrons. Also, states that lie near the continuum 
and have large ionization potentials will be principally populated 
through recombination. when the electron states are close enough to the 
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continuum so the rate of recombination into these states is greater than 
the electron excitation rate. 

As the plasma heats up and highly ionized uranium is present, or 
when a buffer or seed gas is introduced into the plasma the non -Maxwellian 
tail of the electron distribution function can be important in causing 
excitation. To illustrate this consider a 60:40 mixture of uranium to 
helium at a pressure of 1 atm. and a temperature of 5000“K. The 
corresponding electron flux has been calculated by Makowski^^^^ . The 
ground 2j^S and 2^P state will be considered in this calculation, and 
they have densities of 7.82x10^ and 1.26x10^ electrons/cm^, respectively. 
The line temperature [Eq. (V-8)] for the 2^P to ground state transition 
is 4.82x10^ ®K which is much greater than the equilibrium value of 
5000®K. Hence, recombination and the non-Maxwellian tail to the 
electron distribution function can produce high levels of radiation at 
a desired frequency corresponding to line transitions in the plasma. 

C. Compu'cational Method 

The con 5 )utational method developed to do these calculations is 
general. The rate equation formalism can treat both the non-linear 
problem present in this work as well as the "coronal" case. The 
radiative transport problem in cylinder can be used to calculate the 
radiation anywhere on the side of the cylinder. With modifications 
the computer code can treat radiation from any geometric configuration. 
Included here is the radiation from the end of a cylinder which would 
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allow the radiative -energy current from a typical laser configuration 
to be calculated. Also by using an appropriate electron energy 
structure, atomic lifetimes and electron excitation cross-sections, 
other plasmas can be treated. 

D. Future Work 

Due to gaps in knowledge about atomic cross-sections and lifetimes, 
experimental studies should be initiated to measure these. As data 
becomes available it should be incorporated into the uranium model so 
a calculation of exact inversions in the plasma can be made. Studies 
to determine the electron distribution function in uranium mixtures 
should be undertaken so the radiative-energy current can be calculated 
when a buffer or seed gas is present in the plasma. This would allow 
consideration of various radiation processing schemes to be evaluated. 

Also, studies of different plasmas should be undertaken. Impurities 
in a fusion plasma could be investigated as well as the radiation due 
to a neutral gas blanket in a tokamak. The code could be used to 
investigate different lasing schemes both to see if lasing is possible 
and to check the accuracy of models used to explain the mechanism 
driving a particular laser. 
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APPENDIX A. RADIATIVE TRANSPORT CONSIDERATIONS 


A. Review 

In this appendix the expressions for the total radiative-energy 
intensity and the radiative-energy current are derived. Then analy- 
tic solutions for the one-region problem in a sphere and slab are 
obtained. 

The equation for the directed radiative -energy intensity along 
path length s at the origin due to a directed radiative-energy 
intensity at s derived in Chapter II is 


I(v,s) = I(v,o) exp 


I 


Y(v,s')ds' 


+ exp 



s 

Y(v,s')ds* 



f(v.s^) 


ds' 


exp 



YCv,x')dx' 


(A-1) 


where f(V,s) is the volumetric source and Y(V»s) is the attenuation 
coefficient within the voliame. 

When the angular dependence of the gain and the volumetric source 
are known, the total radiative- energy intensity at the origin is just 
the sum of the directed radiative energy flux from all directions, and 
can be written as 


F^ro) = 



d(p 



d(cos0) 


ICV,S[0,(|)]). 


(A-2) 


This quantity is substituted into the rate equations to calculate the 
state densities. 


I 


T 
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To obtain the radistive-energy current in a direction fi' in terms 
o£ the directed radiative-energy intensity at U, the component of the 

/V 

directed radiative -energy intensity in the direction parallel to fi' 
must be obtained; this can be written as 

I(V, s [6.'l>]) = ^ »<{>])• CA-3) 

Now the net radiative-energy current can be obtained by integrating 
Eq. (A- 3) over all angles to give 


- f '■ 

2tt 

1 

; " 

= / d({) 

/ d(cose) 6 • I(v,s[0,(J)]). (A-4) 

' 1 ' 

Jo 

J-1 


This quantity when evaluated at the surface for a given direction 
will give the signal measured by a detector some distance away from the 
plasma. Equation (A- 2) and Eq. (A- 4) can be found in any standard 
reference such as Ref. 5. To obtain the optical power radiated from 
the plasma the current normal to the surface should be integrated over 
all frequencies 


■/ 


J^(o,fi')dv. 


(A-5) 


Simple one region calculations will be carried out for a sphere and 
a slab analytically to illustrate the accuracy of this formulation 

B. One Region Problems 
1 . Sphere 

In this section the total radiative -energy intensity and the 
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radiative-energy current will be calculated in a sphere and a slab as 

illustrations of the use of this form. 

If a spherical plasma of radius a with its center displayed a 
distance or up the z axis (see Fig. A-1) the path length can be 
written as 


s = 2a COS0 . 


(A-6) 


Hence if there is not radiation incident upon the sphere and y and 
f are constant with the plasma the directed intensity for a given 
frequency V can be written as 


f , 2ay COS0 .. 
I(Vr^J=-(e ’ -1). 


(A-7) 


Substituting Eq. (A-7) into Eq. (A-4) with = cos6. The 

radiative-energy current becomes; 

2-n 1 

J^(o,-k) = y* d({) y* d(cos0) COS0 ^ (e^^^ -1), (A-8) 

0 -1 

but f = 0 for COS 0 < 0. Thus 

1 

J,,(o,-lc) = — j d(cos0) COS0 (e^^^ -1), (A-9) 

“ 0 


which upon integration yields 






106 




2tt£ 

Y 


g2aY 

2ay 


(2aY)^ 





2TTf 

Y 


2aY-l 2ay l-2a^y^ 

2 2 2 “ 2 — 

4a^y^ 4a 


When the source is due t« atomic line radiation 


f ANhvgCvj 
4it 


9 


(A- 10) 


(A-11) 


where A is the inverse lifetime, N the emitting state density, v the 
frequency of the light from the atomic transition and g(v) the line 
shape. Substituting this into Eq. (A- 10) the radiative-energy current 
can be written 


8a Y 

This result can also be obtained by considering the geometry in Fig. (A-1) 
and calculating the radiative-energy current going through a surface area 
dot in the x-y plane at the origin. The total radiative- energy current 
passing throusjh da, J^(0-k)da, will be the result of the product of the 
probability of emission of the radiation. 


(2ay-l) 


* Cl-2a\2) 


(A-12) 


prob of emission = AN, 


(A- 13) 
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times the attenuation 


Yr 

attenuation = e' , 


(A- 14) 


times the surface area seen at r of da 


observed surface area 


da cos 0 


4TTr 


2 


(A-15) 


times the probability it is at the right frequency 


line shape = g(v) , 


(A- 16) 


times the energy of the transition 


E = hv , 


(A-17) 


integrated over all space 

= "Mi 

which upon doing the <|) and r integrations yields 

ll ^ ■ 

d(cos0) COS0 , CA-19) 

illustrates the accuracy of the method 


which is just Eq. (A- 9) . This 


2ir 


I 


2a COS0 


d(J> f d(cos0) I COS0 e'^^dr, (A-18) 


Jo 


1 
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developed in Section A. 

For the total radiative- energy intensity, Eq. (A-7) is substituted 
into Eq. (A-2) to obtain: 

2tt 1 





d(cos9) 


f 

a 


j^aaycose 


(A-20) 


which upon integration yields 


F^(o) 



(A-21) 


Another test of the accuracy of these results is to see if the total 
radiative-energy intensity is always greater than or equal to the 
radiative-energy current. To test this in the case of a spherical 
plasma divide Eq. (A-10) by Eq. (A-21) to obtain: 


= -1 




2ya 


(2aY-l)e^^^ h- 1 - 2a V 
e^^^ + 2ay - 1 


(A-22) 


If Y becomes very large this ratio goes to 1; if y a becomes very small 

2 

the ratio becomes -j* 


2. Infinite Slab 

Another simple geometry that can be treated is a slab infinite in 
the x-y plane and with a thickness £ in the z direction (see Figure A-2) 
and having constant physical properties. In this case absorption must 
be considered for analytic results because otherwise infinite intensities 
would arise due to infinite path lengths. Therefore let y = - a then 
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the path length can be written 

s= it/cos0, (A-23) 

substituting Eq. (A-23) into Eq. (A-1) where there is no incident 
intensity, the directed radiative- energy intensity becomes 


Kv.O) = i (A-24) 

vX 

substituting Eq. (A-24) into Eq. (A-2) the total radiative- energy 
intensity becomes 


27T 


F^Co) 


= 1 J d* I 


dCcosB) 


(A-2S) 


which upon integration yields 


2irf 


^ P-E2C«11)] . 


(A-26) 


where E 2 (a^<) is the exponential integral defined by 


E2(ail) 


. cc 

•/ 


•xila 


dx . 


(A-27) 


Relationships and values for the exponential function can be found in 
Ref. 5. 

Substituting Eq. (A-24) into Eq. (A-4) the radiative-energy current 


Ill 


becomes 


2tt 


J^(o,-k) = ^ / d<f> I d(cos0) COS0 (1- 

^0 Jo 


^ -aA/cos0. 

6 J » 


which upon integration yields 




where E^CaA) is the exponential integral defined by 


EjCaA) 


■I 


-aZx 


dx 


To see that the total radiative-energy intensity is greater than 
radiative energy current, take the ratio of Eq. (A-29) to Eq. (A- 
to give 


J^(o.-k) 


j - E^CaH) 
1-E2(a)l) 


but E 2 (ai,) < 2 E 3 (ail) so 


J^(o,-k) 


J - EjCail) 


(A-28) 

(A-29) 

(A- 30) 

the 

26) 

(A-31) 


< 


1 - 2E3(a£) 


(A-32) 
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J^(o,-k) ^ (A-33) 

F^CO) ^ 2 • 

so that the total radiative- energy intensity is always greater than or 
equal to the radiative- energy current. This is because the current 
is just the component of the total radiative intensity normal to the 
surface . 

The accuracy of the preceding analysis depends upon Eq. CA~6) * 

Eq. (A-7), Eq. (A-23), and Eq. (A-24j . These are verified in Ref. 5 
by substituting the formulae found on page 283 into the expression on 
page 22. These results can be used in cases where the energy spectrum 
of the electron flux does not significantly charge over the volume of 
the gas. Then these analytical results can be used instead of a transport 
calculation if the geometry is the same as used in either of the 


examples. 
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APPENDIX B, URANIUM DATA 

In this appendix the optical data for neutral and singly ionized 
uranium is tabulated All state energies are in cm ^ 

(lev = 8066.1 cm" ^). 

The tables are grouped according to the species of atom, the first 

three deal with neutral uranitim and the second three with singly ionized 

uranium. Relative inverse lifetimes for all identified states are found 

in the first table of each group. All relative lifetimes in U° are 

compared to the relative A value of one for the 26226 cm ^ to 0 cm ^ 

(35) 

transition in neutral uranium which has an absolute lifetime^ ^ of 

7- 1 + 

9.785 X 10 sec , while the values for U are normalized to the 

8 - 1 - 1 

2.190 X 10 sec inverse lifetime for the 26191 cm transitions. 

The second table in each group contains the composite state 
relative A coefficients to be used in the calculations, while the 
final table in the group enumerates the sum of the A values that 
give the average A values of the composite state. 




Final 
InitiaiV State 
State 


TABLE B-1. Relative A Coefficients for Neutral Uranlua. 

State energies are in c«-i. 

620 3801 3S6B <276 4iB} 5762 5991 6249 7006- 7104 7326 76«(> 8119 10347 10656 13^ 


13463 

14644 

15632 

15638 

15771 

.00178 

.00642 

.00457 

.00143 

00146 

.00094 

.00281 

.0004 

.00127 

.00101 

.00203 


.00047 

.00075 

..00108 

.00266 

.00096 

.00159 

.00213 

.0C014 

.00157 

.00296 

.00255 

.00316 

.00137 

.00239 

.00348 

.00166 

.00237 

.00356 

.00103 


.00047 

.00370 

.00319 

.00152 


16122 


.00554 




.0039 

.00744 

.00995 









16506 

.00097 

.00109 

.00091 


.00096 




.00037 

.00040 




.00034 



16900 

.01511 


.00082 


.0006 




.03065 

.27026 


.03290 


.03578 

.01449 

.00731 

16930 

.00054 

.00166 



.0006 

.00020 

.00101 


.00106 

.00113 







17070 

.00221 

.00549 

.00341 


.00057 


. 00095 



.00374 


.00386 


.00421 



17262 

.00169 

.00367 



.00098 



.00192 

.00028 

.00230 


.00153 


.00100 



17468 


.00105 




.00004 

.00062 

.00094 



.00139 






17894 


.00013 





.00019 

.00319 









17908 

.00011 

.00024 



.00007 

.000.38 

.00032 

.00338 

.00039 








17969 




.00035 


.00108 


.00365 









18754 

.00053 

.00105 

.0011 




.00058 


.00089 

.00164 







18406 

.00232 

.00457 



.00205 

.00079 


.00036 

.00376 

.00215 







18531 




.00046 


,00243 





.00111 






18759 

.0041 

.00404 

.00089 


.00044 


.00258 


.00040 




.00468 

.00489 



16933 

.0073 

.00168 



.00161 

.00018 

.00033 

.00144 

.00331 

.00103 







19126 


.00084 


.00044 


.00009 


.00340 









19192 


,00112 




.00019 

.00089 

.00143 



.00089 






19472 

.00014 

.00002 



.00024 

.00016 


.00003 


.00035 







19489 



.00289 










.00089 

.00754 

.00164 

.001 » .00068 



.00123 


.00117 


.00065 

.00007 

.00036 



.00035 






19648 

.00259 


.000329 


.00075 





.00042. 


.00041 



.00053 

.00113 

19783 

.00179 

.00112 

.00060 




.00040 


.00005 

.00088 


.00010 


.00086 



19885 

.02285 


.00088 


.00037 




.00055 

.00064 


.00041 

.00068 

.0004 


.00179 

20114 

.00172 

.00081 



.00043 

.00040 

.00007 

.00020 

.00062 

.00069 







20219 

.0007? 

.CGC31 

.00095 


.00063 


.00053 



.00038 


.onsnn 


.00107 



20312 


.00117 



.00875 




.00169 

.00055 







20392 




.00056 


.00016 


.00005 









20420 

.00057 

.00055 

.00026 


.90046 

- 

.00015 



.00024 


.00053 


.00047 



20464 

.00116 


.00144 


.00192 







.00056 


.00065 

.00093 

.00060 




.00657 










.noDsn 

.00034 

.00477 

.00493 

20569 


.00309 


.00182 


.00419 

.00010 










20621 





.00091 


.00019 


.00033 

.00035 







20662 


.00023 

.00136 


.00658 


.00012 





.00009 


.00117 



20766 

.00015 


.00082 


.00111 




.00012 

.00010 


.00065 

.00116 

.00042 

.00130 

.00135 

20852 

.00012 

.00029 



.00179 

.00043 

.00050 


.00045 








20943 

.00085 

.00050 



.00144 


.00055 



.00004 


.00041 


.00005 
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TABLX B-1. (cciiCiiiucd) 


■ Final 
State 


Initial' 

State 


21062 

21079 

2 1265 

21SES 

21637 

21766 

21768 

21941 

22038 

220 S 6 

22368 

22378 

22464 

22583 

22633 


23187 


24757 


620 3 S 01 


3868 


4276 


4453 5762 


59 P 1 


6249 


7006 


7104 


7326 7646 5119 10347 10686 13128 


.00065 


.00018 

.00533 

.03888 

.08329 


.00036 

.0073 

.00591 

.02394 


.00313 

.06233 


.001 

.00165 .00220 

.00435 

.00474 

.00085 .00048 


.00019 

.00021 

.00293 

.02658 .00075 

.00115 


.00009 


.00271 

.00318 

.01617 .00742 

.00026 .00079 


.00039 


.00213 


.00033 

.00247 


.00036 


.00019 

.00533 


.00002 

.00136 .00883 .00114 


.00016 

.00108 


.00034 

.00023 


.00370 


.00157 

.00107 


.00088 

.00105 
.00277 .00112 

.00138 


.00023 .00000 

.00009 


.00007 


.00005 .00064 

.00093 


.00081 


.00888 .00158 


.01106 

.00004 

.00008 

.00860 

.00009 

.00019 .00064 .00179 

.00072 

.00974 

.00165 

.00307 

.00284 

.00146 

.00002 

.00877 

.00053 

.00458 

v 00416 

.00139 

.00024 

.00067 


.00388 .00022 

.00089 .00010 

.00249 .00060 

.02299 

.00108 . 00091 


.00037 

.00046 

.00016 

.00012 


.00093 .00050 .00026 .00034 
.00104 .00033 

.00012 .00054 .00029 


22754 

.07285 

.01306 

.00461 

.00097 

22790 



.00109 


22862 

.01558 

.00986 

.00782 

.00252 

22919 

.1011 


.03922 

.00921 

23058 

.02942 


.02611 

.00009 


.00016 


.01169 


.00097 


.00023 

.00165 

.00046 

.00128 

.01625 



.09397 

.01525 

.00045 

.04477 

.03318 

.02432 

.02387 

.00878 


.00012 

.00006 

.00187 


.01228 


23197- 

.00067 


.00069 

.00149 

23325 

.00012 

.00001 


.00050 

23430 


.00048 



23433 

.01433 

.01224 


.00801 

23544 

.06995 


.00442 

.00009 

23572 

.04749 

.1019 

.03412 

.00064 

23649 

.01587 


.00024 

.00108 

24067 

.2721 


.00333 

.00160 

24186 

.06465 


.00926 

.01841 

24334 

.0224 


.00533 

.00222 

24433 

.03664 

.03932 

.00424 

.01065 

24560 

,00058 


.00228 

.00586 

24581 



.00565 


24671 

.00292 

.1167 

.00674 

.01580 


.02127 


.02044 


25319 

.1254 

.1467 


.02391 

25349 

.5169 

.87145 

,01836 

.0383.5 

25389 

25463 

.1418 

.68751 

.00013 

.00449 

.02579 


.00108 

.00102 

,00087 

.02227 

.00171 .00081 

.00023 .00084 

.00145 .00147 

.01441 .02160 

.00138 

,00079 

.00372 

.00038 

.02695 

.00007 

.00053 


.00013 .00058 



.00035 


.00499 

.00023 


.00046 

.00008 . OOOC 9 

.00017 




.00084 

.00103 


.00166 






.00107 

.00131 


.00056 

.00026 

.00013 

.00033 



.00469 



.00463 

.00375 .00074 





.00097 

.00065 


.09063 

.26304 


.01270 



.00798 



.00456 

.00173 

.00139 



.01741 

.05735 

.00172 


.02325 

.00093 



.00029 


.00159 

.00228 


.00270 

.00953 

.00104 

.00039 






.00277 

.00230 .01031 

.00023 

.00076 


. 10269 

.05793 

.06474 


.00083 





.00503 .02329 



.01877 






.00902 

.03994 

.04775 







.01568 

.31357 

.03191 


.06943 

.00287 






. 


.00104 

.00231 .00065 


.00039 


.01348 

.05296 

.00448 


.01315 

.03151 
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\ Final 
\ State 
InitialST 
State N 

0 

620 

380! 

3868 

TABLE 8-1. (continiici!) 

4276 4453 5762 5991 6249 7006 7104 

7326 

704f5 

8119 

10547 

10686 

13128 

2S627 





. 


















,00057 

.00020 

.00058 

2S672 

.2240 


.00977 


.00660 

.07437 


.04428 

.07859 

.00660 



2S826 

.7112 

.03737 

.04799 


.00039 .00671 

.04288 .0X817 



.09681 




25918 



.00043 




.00100 

.0022 

.00112 

.00035 

.00019 

.00049 

25938 

.3058 

.05444 

.04297 


.01739 .05065 


.07570 


.02050 




26104 

.06236 

.02164 

.00878 


.00592 .00401 

.00305 .00725 

.02805 






26226 

1.0000 

.51332 

.08091 



.14962 







26313 



.00859 




.11694 

.00519 

.00884 

.01634 

.00712 

.00321 

25454 



.00549 




.00554 

.01588 

.00916 

.01170 


.00617 

26550 

.1549 

.09340 

.00534 


.00386 

.00845 .01239 







26608 

.U0032 


.0123 


.00205 

.00619 .00359 


..00172 

.00326 

.00568 

.00123 


26792 

.05759 

.02402 

.00212 


.00164 .00039 

.02117 .01629 







26921 

.2209 

.22611 

.05759 


.02814 .01632 .04062 .04503 

.01442 .03674 







26979 



.03703 




.02400 

.00012 

.00088 

.00129 


.00013 

27524 

.00389 


.00064 


.00421 

. .00607 .00155 

.00039 

.00410 


.00161 

.00824 


27478 ■ 



.09889 




.00093 


.00672 

.00320 


.00052 

27499 

.06136 

.06726 



.00094 .01479 

.00446 . 00249 



.00162 




27616 

.1315 

.1471 

.00636 


.00592 .00381 

.00151 .00306 







27682 

.03086 

.01477 



.00504 .00639 .00743 .00065 

.00195 .00238 







27744 

.06267 

.09834 

.02047 


.02510 

.00091 







27791 

.02i.27 

.1266 



.01759 .00089 

.00593 .00012 







27818 



.07442 




.00575 

.00571 

.00386 

.00041 

.00228 

.00176 

27887 

1.407 


.00576 



.05094 .00341 



.00584 




27941 

.1646 

.5192 

.04466 


.05112 

.01433 .00679 







27966 

.1436 


.05762 



.01140 .01029 

.00849 


.01932 

.03075 



28053 

.1806 

.25229 

.07585 


.00297 

.00286 . 00167 

.00016 


.02542 




28099 

.1827 

.09581 

.00058 



.00170 .02994 







28119 

.3214 

.01766 

.02190 


.00354 

.01848 .00145 

.00466 


.00463 . 

.01766 



28153 



,09099 


.05592 

.00495 .00115 

.01686 



.06382 

.40117 


28286 

.1556 


.01034 


.02473 

.03093 .02855 

.00960 

.00182 


■00116 

.00605 

-00697 

28444 

.8097 

. 1705 



.00410 .00399 

.00557 .00530 







28454 


.0641 


.06036 

.00068 

.00176 .00107 







28470 

.1254 

.1073 . 

00347 


.00110 .00109 

.00323 .00086 







28504 • 

.4047 

.13387 



.00928 .00514 .00140 .00776 

.00732 .00229 







28543 

.09162 

.0874 . 

00784 . 

.04626 

.02007 .00234 .00384 

.00275 .01133 







28563 

,3338 

.08770 



.-00682 , 00058 

.00098 







28650 

.8529 1.0662 



.00364 .00184 .00165 .00063 

.01063 







28799 

.00441 

• 

04739 


.05808 

.00814 

.00310 . 

.00281 . 

01170 . 

0S9S9 



28875 

.1806 

.5939 . 

03963 


.08440 .00785 

.02771 .00657 

.00375 






29126 

.00516 . 

.1365 . 

00472 


.02408 .00022 

.01652 .00121 
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X Finil 
NStale 
InltlalX 

0 

620 

3801 

3868 

4276 

4453 

TABLE B 
5762 

-1. (ciHilinucd) 
5991 6249 

7006 

7104 7326 

7646 

8119 

10347 106S6 

15128 

State ^ 















29233 


,1811 




.07583 

.02535 

.00410 

.00396 






29414 



.02836 



.1670 




.05443 

.01593 

.0.8160 

.08505 .14706 


29559 

,023’3 

.21543 

.00216 





.05215 


.00738 


.00567 



29791 


.08795 

.1212 


.02962 

.00462 

.00627 

.00414 

.03025 

.01573 

.00090 

.00746 



29S38 

1.66273 


.6133 


.38725 



.02055 

.01982 


.10401 


.06659 


29S66 

.02764 

.04577 

.02758 


.01772 


.00271 

.02256 




.00182 



•• 29958 

.25831 


.1262 


.00148 



.02059 



.00955 

.01512 



30222 






.1148 









30451 

.00671 

.06249 

.3439 


.06347 


.04446 

.02733 

.00339 

.00525 

.00817 




30499 


.1075 


.05145 


.00853 

.00033 

.00353 







30587 

.00293 


,1397 


,00729 



.00140 

.00934 

.02550 


.01263 

.01838 


30637 

.1086 


.04967 


.12702 


.03504 


.01242 

.02614 


.00620 



30876 

.02682 


.04126 


.2823 


.04877 

.01980 

.00347 



.00448 



30986 



.1362 







.03122 

.01586 

.02184 

.02193 


30993 

.00883 

.01333 



.00490 

.2474 

.04124 

.00512 .00667 







31221 


.03787 


.07579 


.01409 

.1417 








31270 











.06007 


.04586 .03152 

.07822 

■ 31276 

.00900 


.04080 


.1473 



.03397 


.04901 

.00154 

.09013 

.00388 


31279 



.4689 







.22424 

.00584 

.31284 

.38838 


' 31301 



.07132 


.03096 



.00124 

.00031 

.00493 

.00535 


.002406 


31445 



.1324 







.00283 

.00212 

.02604 

.070857 


31552 



.06278 


.04318 


.01022 

.00097 

.01352 

.00300 





31946’ 

.09091 




.10496 

.36627 

.71612 

.06009 .12976 







32017 

.00370 


.04468 


.05218 



.00870 

.01863 

.01542 

.00053 

.04923 

.00732 .00252 


32098 

.2208 


.12882 


.106$ 



.05004 

.01387 

.02678 . 

.00567 

.05915 

.00431 .01125 


32318 

.09405 




.06754 

.05384 


.03895 .02194 







32492 

.00812 

.00429 



.1531 


.02920 

.2011 







33412 

.00407 

.4523 



.08902 


.00960 

.00224 

.01668 

.01341 


.01824 



S3S71 

.03625 


.17039 


.11180 



.03153 

.1866 


.02961 

.00224 

.03306 .0121$ 


33640 

.03189 

.3997 

.18937 





.05950 

.00641 

.08134 


.1182 
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Initial' 


Final 

State 


620 


3801* 


4276 


44S3 


TABLE B-2. (continued) 


5762 


24450 

.10158 

.15602 

.02024 

.04156 

24558 

.02298 

.02127 

.03370 

.00808 

25451 

1.0081 

1.1057 

.03262 

.09464 

25645 



.00061 


25997 

2.017 

.60514 

.17187 

.01778 

26209 

.06236 

.02164 

.01737 

.00592 

26531 

.00032 


.01779 

.00205 

26754 

.43339 

.34352 

.06505 

.02978 

27392 

.09811 

.08203 

57204 

.13656 

.15887 

.01018 

.04861 

27914 

1.5716 

.5192 

.10228 

.05117 

2E192 

1.6500 

.53627 

.19967 

.08829 

28556 

.95992 

.48037 

.16532 

.09535 


.00029 


.12010 

.00503 

.02916 


.05065 


5991 

6249 

7006* 

7326 

7646 

8119 

U347 

10686 

15128 


,09624 

.06710 

.11471 


.27023 


.01270 


.02329 

.00957 

.02106 

.01003 

.00230 

.03157 

.00266 

.00039 

.00076 

.00903 

.48034 

.08414 

.08258 

.04428 

.11297 

.00860 






.00204 

.00253 

.00177 

.00092 

.00038 

.00146 

.00671 

.04288 

. 16779 

.07370 


.11731 



- 


28763 

29425 


.01632 

.00639 

.00089 


.00401 

.04488 .04503 
.02222 .00065 
.00381 


1.0335 

.02829 


1.6601 

.61098 


.03963 

.15643 


.08804 

.05370 


.00582 

.00184 

.25145 


.00696 

.00541 .01160 
.00951 .00063 
.03185 


00305 

.00725 

. 14499 

.00519 

.00884 

.01634 

.00712 

.00321 

00619 

.00359 

.00554 

.01760 

.01242 

.01758 

.00123 

.00617 

04403 

01249 

.06542 

.00641 

.02532 

.00422 

.00922 

.00610 

.00824 

.00066 

01976 

.01347 

.01425 

.00571 

.02319 

.03115 

.00228 

.00176 


.06527 

.06450 

.01020 

.06807 

.03128 

.00182 

.00584 

.03004 

.08263 

.40722 

.00697 

.01758 

.03615 

.00310 

.00281 

.01170 

.0S505 



.03834 

.07891 

.00657 

.03541 

.00375 

.07754 

.01683 

.09473 

.14706 



29838 

1.6627 


.6133 

.38725 




.02055 

30015 

.28594 

.04577 

.15378 

.01920 

.1148 

.00271 


.04316 

30718 

. 15390 

.18332 

.76218 

.48502 

.25593 

.16985 

.00865 

.05320 

31344 

.00900 

.03737 

.38309 

.22144 

.01409 

.15192 


.03619 

31613 

.09091 


.4689 

.10496 

.36627 

.71612 

.06009 

.12976 


32231 

33541 


.23675 

.07221 


.00429 

.852 


.17351 

.35976 


.37932 

.20082 


.05384 


.02920 

.00960 


.03895 


.01982 

.02022 

.01383 


.09111 

.05976 

.22723 


.28178 

.09326 


.03271 

.20969 


.10401 .06659 

.00955 .01693 

. 02404 . 04515 . 04030 

.06908 .11617 .07300 .03152 

.00584 .31284 .38838 


.04220 

.09476 


.00620 

.02961 


.10838 

.13667 


.078222 


.01163 

.03306 


.01377. 

.01215 
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\ Final 
\ State 

InitialV ® 

State ' 

TABLE B.3. 

620 3801» 4276 4453 

U* Kelative Coefficients 

State energies arc in cia 

5762 5991 6249 7006* 

7326 

7646 

8119 

10347 

106S6 

13128 

<k> 

*108 

13463 

.00178 

.00094 .00047 

.00108 

.00096 

.00014 








.00526 

14644 

.00642 

.00281 




.00157 

.00165 


.00047 




.01265 

1S778 

.00249 

.00240 .00152 .00075 

.00152 

.00372 

.00995 

.00289 .00241 

.002.32 


.00344 

.00152 



.01808 

16S06 

.00097 

.00109 .00091 .00096 




.00037 .00041 



.00034 




.09493 

16900 

■OlSll 

.00082 .0006 




.03065 .27026 

.03291 


.03578 

.01449 

.00731 


.35913 

17207 

.00148 

.00297 .00341 .00072 

.00012 

.00086 

.00144 

.00067 .00208 

.00280 


.00261 




.01176 

17924 

.00011 

.00019 .00035 .00007 

.00073 

.00025 

.00041 

.00039 







.00146 

18397 

.00143 

.00208 .00078 .00205 

.00161 

.00058 

.00037 

.00232 .00164 







.00842 

18846 

.0057 

.00286 .00089 .00103 

.00018 

.00146 

.00144 

.00186 .00103 


.00468 

.00489 




.01904 

19263 

.00014 

.00067 .00044 .00024 

.00015 

.00089 

.00063 

.00057 







.00233 

19618 

.00219 

.00117 .00125 .00075 

.00066 

.00024 

.00036 

.00046 .00055 

.00025 

.00089 

.00081 

.00109 

.00117 

.00068 

.00587 

1988S 

.02285 

.00088 .00037 




.00055 .00064 

.00041 

.00068 

.00040 


.00179 


.02795 

2021S 

.00124 

.00077 .00095 .00.327 

.00040 

.00030 

. 00020 

.00117 .00054 

.00504 


.00107 




.00874 

20406 

.00057 

.00055 .00041 .00046 

.00016 

.00015 

.00005 

' .00024 

.00053 


.0004? 




.00194 

205S6 

.00116 

.00166 .00280 .00425 

.00419 

.00011 



.00038 

.00050 

.00072 

.00285 

.00276 


.01057 

208S2 

.00037 

.00040 .00075 .00131 

.00072 

.00042 

.00019 

.00030 .00014 

.00053 

.00116 

.00024 

.00130 

.00135 


.00360 

21172 

.00425 

.00383 .00313 .00300 

.00220 

.00157 


.00309. 

.00039 


.00213 




.01925 

21611 

.06109 

.01493 .06233 .00279 

.00048 

.02658 

.00095 

.00822 .00410 







.13382 

21878 

.00140 

.00018 .00182 .00028 

.00132 

.00182 

.00108 

.00016 .00041 

.00007 


.00005 

.00788 



.00346 

22S1S 

.01717 

.01005 .00473 .00239 


.00384 

.00138 

.00675 .00406 

.00030 

.00098 

.00071 

.00026 

.00034 


.04498 

22600 

.00003 

.00030 .00600 .00095 

.00019 

.00064 

.00018 

.00099 .00051 • 

.00088 

.00029 

.00030 

.00043 



.00459 

22989 

.06526 

.03267 .00465 




.06358 .01978 

.01216 

.02678 

.00029 


.01226 


.22618 

2328S 

.00039 

.00024 .00068 .00100 

.00099 

.00113 

.00100 

.00108 .00118 


.00013 

.00058 




.00514 

23691 

.11131 

.05707 .01872 .00342 

.02228 

.01441 

.02160 

.00276 .01399 

.00315 

.00375 

.00074 




.19281 


s 




« 
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Initial' 

State 


3801* 4276 44SI 5762 


TABLE B-3. (continued) 

5991 6249 7006* 7326 7646 8119 10347 106S6 


-04251 .00233 .00038 

.03386 .07801 .00675 .01496 

.01149 .02127 .00843 .00404 

.25202 .36858 .01087 .02366 

.00031 

.67233 ,20171 .05729 .00889 

.06236 .02164 .00868 .00592 

•00032 .00889 .00205 

.14446 .11451 .02168 .01489 

.03270 ,04101 .04552 .00339 

.09076 .12401 .03972 .01620 

.7858 .2596 .05114 .05117 


.00038 .00035 

.01496 .06005 

.00404 .00029 .00303 


.00303 .'00077 
.03208 .02237 
.00479 .01053 
.12021 .02805 


.00030 . 00008 
.03624 

.00334 .00230 

.04129 . 04428 

.00102 .00126 
.07370 


28192 

.33 

.13407 

.03993 

.02207 


28556 

- .19198 

.09607 

.03306 

.01905 

.00291 

28763 

.51675 

.83005 

.03963 

.04402 

.00184 

29425 

.01414 

.15275 

.03911 

.02685 

.08382 

29838 

1.6627 


.6133 

.38725 


30015 

.14297 

.04577 

.07689 

.00960 

.1148 

30718 

.03078 

.06111 

. 12703 

,097(W' 

.12797 

31344 

.00900 

.03787 

.07662 

.07381 

.01409 

31613 

.09091 


.4689 

.10496 

.36627 

32231 

.05919 

.00429 

.01675 

.09483 

.05384 

33541 

.02407 

.426 

.17988 

.10041 



.05065 .00671 .04288 .08390 .07370 

.00401 .00305 .00725 .07249 ,00519 

.00619 .00359 .00554 .00880 

.01632 .01496 .04503 .01468 .02181 

.00639 .01111 .00065 .00416 .00214 .00844 .00211 

.00089 .00381 .00494 .00449 .00712 .00571 

.03214 .00510 

.00348 .01075 ,01134 .00782 .00182 

.00291 .00135 .00580 .00293 .00904 .00310 .00281 

.00184 ,00475 .00063 .01917 .00657 .00 375 

.08382 .01062 .01973 .01180 .02585 .00841 

.02055 .01982 .10401 

.1148 .00271 .02158 . .00955 

.12797 .03397 .00432 .01380 .00505 .02278 .01202 

,01409 . 07596 .01206 . 00692 . 01494 . 01727 

36627 .71612 .06009 .12976 .11362 .00584 

05384 .02920 .03895 .07045 .01635 .02110 .00310 


.00017 

.00015 

.00C55 


.09008 


.01270 


.010522 

.00039 

.00039 

.00076 

.03766 

.00860 



.00089 

.00046 

.00016 

.00049 

.03867 




.00884 

.01634 

.00712 

.00321 

.00621 

.00879 

.00123 

.00617 

.00307 

.00203 

.00824 

.00033 

.01 159 

.01558 



.00584 




.01502 

.02754 

.20361 

.00697 

.01170 

.05959 




.00841 .05158 .08505 

.10401 .06659 

.00955 .00847 
.01202 .01129 

.01727 .05809 .01825 
.00584 .31284 .38838 

.00510 .05419 .00582 

,P2961 .04556 . 03506 


.04975 
.32739 
.05360 
.07565 
■00517 
1 .06680 
. 15054 
.04426 
. 35467 
.33912 
.20720 
1.10981 
.51756 
.29584 
1.40993 
.32647 
2.81238 
.22565 
.32077 
.21122 
1.40528 
.34553 
.68610 
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\ Final 
\ State 
InitialX 
Stftte ^ 

0 

289 

915 

TABLE B-4. 
1749 2295 

Relative A Coefficients for Sinclcly Ionized Uranium. 
State energies arc in cm'^. 

4421 4585 4706 5260 5402 5526 5667 

5791 

6283 

6445 

8394 

1S200 

.00605 

.00162 

.00192 


.00251 












20572 

.00451 


.00183 





.00139 


.0020 


.00124 





20962 

.00499 


.00234 





.00623 


.00024 


.00134 





21021 

.00674 


.00237 





.00674 


.00017 


.00018 



.00056 


21320 

.01239 


.00981 







.00064 


.00832 



.00016 


21711 ~ 


.01020 


.00654 



.00169 


.00354 




.00026 




22101 

.00032 

.00038 



.00100 

.00077 






.00027 



.00048 


2216S 

.00555 

.01034 

.00159 

.00120 

.00092 

.00491 

.00944 



.00260 


.00047 

.00403 




22430 

.00270 

.01985 

.01706 



.00870 




.00099 


.01118 



.00036 


22642 

.01245 

.05108 

.00399 


.00430 

.00293 






.02838 

.00044 


.07049 


2286B 


.00981 

.00316 









.01258 



.00072 


22917 

.01864 

.00294 

.04905 

.01134 

.00459 

.00570 

.01905 








.01015 


23241 

.01297 

.00646 

.00577 

.00126 

.00661 

.00175 

.00055 



.00781 



.01513 

.00817 

.00056 


23315 

.02201 

.06837 

.01292 


.00027 





.000728 


.01123 


.01761 



23554 

. 13167 

.01595 

.01364 

.00407 

.00713 

.00106 





.00651 


.00047 

.00070 

.00067 


23636 


.02218 


.00726 

.00582 


.00211 


.00042 





.00459 



23778 

.02137 


.01577 

.01503 

.00329 

.00257 







.00082 


.01253 


23912 

.00361 

.00791 

.01494 



.00377 




.00139 


.01333 



.00128 


24010 

.04624 

.00972 

.00492 

.00908 

.02133 

.00733 

.00008 






.00130 


.00131 


24153 

.02323 

.03693 

.00253 

.00032 

.01236 

.00305 







.00135 

.00095 



24160 ~ 


.02354 


.01582 

.00742 

.03423 



.00357 





.00790 


.00470 

24288 

.11178 

.03470 

.00092 

.01049 

.01236 


.00549 




.00258 



.00066 

.00080 


• 24506 

.01497 

.00269 

.00897 


.00022 

.00410 




.00082 


.00077 



.00100 


24342 

.02332 


.01508 





.01216 


.00163 





.00621 


24453 

.03881 

.00019 

.01136 


.00428 

.00452 




.00084 


.00069 



.00103 


24608 

.13103 

.07989 

.00986 

.01050 

.00183 


.00379 






.00070 

.01251 

.00150 


24684 

.24648 

.06505 

.00038 


.00586 

.00176 




.00278 


.00209 

.00577 


.00262 


25164 


.04737 


.01804 


.00315 

.00530 


.00135 


.00072 


.00398 

.00022 


.00111 

25357 

.03397 

.15349 


. 16678 

.02199 

.06828 











25714 

.00082 

.03701 

.02686 

.17207 

.01331 

.00403 

.00864 


.00088 

.00168 



.00518 

.00115 



25968 

.07213 


.07651 


.01092 



.01841 

.00640 

.03636 


.04521 


.00043 



25986 

.44247 

.04236 


.05519 



.00738 

.00357 

,00015 


.00184 


.00173 

.00042 



26191 


1.0000 


.3096 

.00795 

.00152 

.01842 


.00117 


.00149 


.00501 

.00576 



26415 


.16154 


.02658 

.01963 

.00493 

.01087 




.00639 


.00296 



.00199 

26717 


.44334 


.05729 

.00294 


.01821 







.01059 

.00868 


26887 ~ 







.00132 


.00291 





.00279 


.00093 

27126 

.00018 

.00086 

.08400 


.01278 




- 

.00518 


.00105 

.00525 


.00082 


• 27290 

.00017 


.01953 





.00614 


.00033 


.00217 



.00055 


27357 

.03931 

.03156 

.08898 

.01328 

.03157 

.00677 

.01880 






.00628 

.00287 

.00678 


27390 



.05767 

.02680 







.01404 
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table B-4. (continued) 


\ Final 
. , \state 
Inltlai\ 

0 

289 

SIS 

1749 

2295 4421 

458S 

4706 

5260 

5402 

5527 

5667 

5791 

6263 

6445 

8394 

State ' 

27499 

27930 

28154 

28308 

.02684 

.02350 

.23431 

.05460 

.01491 . 

.12754 

.00197 

.14669 

.62061 

.06783 

.09629 . 

.02680 . 
.00455 , 

.01717 

,04205 .00822 

.11957 .02070 

.38475 .01626 

.08404 .00269 

.01808 

.00665 

.01774 

.11714 

.01704 


.00069 


.05974 . 

.00571 

.00205 

00597 

.00076 

.00711 

.00257 

.00022 

.00459 

.00310 

.00221 

- 

.00545 

28578 

28758 

.02480 

.01526 

.03171 

.05868 


.00381 .00150 

.01268 




.00475 

.00532 



.00038 

.00575 



29207 

29477 

29828 

.00922 

.02297 

.02877 

.04684 

,06234 

.05868 

.01083 

.05737 

.02773 

.02874 

.06076 .02578 
.00258 .02123 
,00005 .00272 

.01704 

.04357 

.03290 

.00098 

.08192 


.02743 

.00178 


.00348 

.00143 

.00570 

.00116 

.00796 


.00220 

29932 

29956 

30061 

30C86 

30240 

30264 

30342 

.02834 

.04745 

.00916 

.36416 

.06407 

.08084 

.02219 

.12304 

.00718 

.00571 

.00882 

.02672 

.05941 

.07868 

.04922 

.08039 

.00022 

.06908 

.01514 

.00628 . 03644 ’ 

.00352 .03928 

.00907 .00006 

.02464 

,00988 

.03499 

.00065 

.01883 



.00267 

.00647 

.00501 

.00324 

.COOS? 


.00009 

.0273C 

.00237 

.00232 

.00661$ 

.01717 

.00013 

.00028 





.11996 

.00229 

.00840 

.02312 .01768 

.05418 

.201 


.19264 


.00091 






30469 

305S0 

30860 

30863 

30942 

31084 

31219 

32211 

.00143 

.08805 

.16207 

.02086 

.02060 

.04607 

.00165 

.0683? 

.00335 

.00276 

.01735 

.01743 

.04642 

.01422 

.00927 


.00895 

.00023 

.00760 

.07602 

.00030 

.00855 

.00053 

,00271 

.08579 

,27626 


.1523 

.03077 

.02409 

.03438 .02186 

.00331 
.1372 
.49839 

.1307 

.00205 

.00328 

.00118 

.00110 


.12098 


.00184 
.00573 
.00959 
. 16879 
.29286 


.00049 

.01934 

.10113 

.15078 

.04491 
.00796 
.01967 
. 00275 
.00229 

.00027 

,00609 

.00355 

.02157 

.06109 

32535 

33104 

3347S 

33749 

34049 

34439 

34708 

.01140 

.12788 

.07034 

.13205 

.00906 


.09897 

.07255 

.01000 

.14006 

.05740 .01552 

.06310 

.01223 

.00731 .02075 

.00026 .00097 

.02651 

.00458 

.02040 

.00077 


. ,05418 
.04953 

.06601 


.00755 

.01442 

.00619 


.01029 

.00250 

.00432 

.00339 

.05750 

.01076 


.00207 

.00655 

.09405 


.24164 


.35175 

.07210 

.02863 


.09477 





.03136 


•0SS40 

34886 

34924 

1.0906 

1.2651 

20441 

.23678 

.20837 

.00025 

.06090 . 

.00986 

.00019 


.00305 


.00038 


.0375 

.00039 


.00001 

3S6ES 

55790 

36789 

37308 


.497B1 

.23813 

.19199 

.08767 

.00692 

.41889 

.01908 


.00025 


.01009 


.02335 

- 



.41373 


.22663 

.14759 

.00779 


.06888 

.28886 


.00481 

.24525 






391C9 

S9S0S 

1.3096 

.05016 


,41466 

.26978 

.11014 


.16184 


.1187 



.00112 

.03102 
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table B-S. Radiative Relative A Coefficients for U . 
State energies are in c» . 


Final 
, State 


Initial' 

State 


289 


91S 


1749 


2295 


24183 

24367 

24818 

2S680 

26091 

26303 


18200 

.00605 

.00162 

.00192 

20S51 

.00483 


.00218 

21516 

.01259 

.01020 

.00981 

22101 

.00032 

.00038 


22298 

.00413 

.01510 

.00933 


.00251 


.00654 


.00120 


.00100 

.00092 


22780 

.01545 

.02701 

.02652 

.01134 

22868 


.00981 

.00316 


23570 

.05555 

.03026 

.01077 

.00276 

23778 

.01249 

.01505 

.01535 

.01114 


.06042 

.02737 

.18876 

.03564 

.42247 


.02622 

.00144 

.06410 

.09523 

.52118 

.30244 


.00279 

.01180 

.00512 

.05168 


.00445 

.00467 

.00456 

.01782 


4421 

4585* 

5330* 

5596* 

5791 

6283 

6445 

8394 

<A> 

xlO* 









.02652 


.00479 

.0002 

.00050 


• 

.00056 


.03006 


.00169 

.00209 

.00832 

.00026 


.00016 


.05864 

.00077 



.00027 



.00048 


.00702 

.00680 

.00944 

.00180 

.00583 

.00403 


.00036 


. 11160 

.00432 

.01905 

.02838 

.00044 



.00054 


.32361 



.01258 




.07049 


.08164 

.00141 

.00055 

.00427 

. 00887 

.00780 

.00885 

.00562 


.27891 

.00317 

.00211 

.00091 

.01335 

.00082 

.00459 

.00691 


.11680 

.01267 

.00278 

.00357 

.00258 

.00132 

.00317 

.00106 

.00470 

.24769 


.01427 

.16942 

.18240 

.04194 


.00230 

.00431 

.01216 

.00110 

.00073 



.00241 


.00384 

.00245 

.00455 

.00206 

.00140 

.00348 

.00637 

.00206 

.00111 

.01541 

.03616 

.01352 

.01133 

.04521 

.00318 

.00079 



.00795 

.00152 

.00979 

.00066 

.00167 

.00337 

.00309 



.01128 

.00493 

.01454 


.00639 

.00296 

.01059 

.00868 

.00199 


.11591 

.48394 

.71728 

2.06537 

.84973 

.04035 

.39334 

2.06903 

.39794 

.39305 


30291 

30627 

3090! 

31152 

32373 

33290 

33749 

34471 

34797 

35738 

37049 

39309 


27088 

.00017 


.01953 




.00373 

.00162 

.00217 

27357 

.02211 

,01577 

.05815 

.04546 

.02879 

.00750 

.01273 

.00518 

.01919 

28042 

.12891 

.12754 

.38365 

.01567 

.25216 

.01848 

.06744 

.00069 

.00571 

28614 

.03970 

.02348 

.05570 

.05419 

.02392 

.00743 

.01486 


.00340 

7Q<;n4 

.01609 

.04599 

.03476 

.03795 

.02113 

.01641 

.01704 

.00098 

'.OllSl 

30061 

.02832 

.05570 

.02028 

.05213 

.00490 

.02887 

.C2440 

.04224 

.00412 


.00517 

.00484 

.00246 


.00279 

.00181 

.00384 

.00130 

.01261 


.00055 

.00261 


.00545 


.00220 


.36416 

.00143 


. 12304 
.05754 
. 16207 


.02672 

.02073 

.04607 


.09452 

.00250 

.11033 


.00907 

.03129 

.02312 


.00006 

.00229 

.01768 


.10992 

.00276 

.05516 


.19264 

.04642 

.06760 


.00208 

.00927 

.00184 

.00776 


.01370 


.00459 


.01292 


.00481 

.02620 

.00760 

.06047 

.01382 


.00013 

.00028 .00631 
.00030 .00053 
.00855 .00609 
.000272 .00355 


.40855 

1.30065 

.22189 

1.94307 

.61006 






.31780 


.00114 


.23083 

.12591 

.00152 

.04133 

1.57597 


.09911 


.08576 

.06025 

.01552 

.02651 

.05186 

.01099 

.01029 

.03044 

.00136 

.6/540 

.01140 

.13205 


.01000 

.01223 


.00(58 



.00250 



- .37839 

.09616 

.17186 


. 14006 

.02284- 

.01086 

.00712 

.03453 

.03282 

.02091 

.00558 

.04703 

.69497 

1.0906 

.75337 

.23678 

.27006 

.07210 

.00025 

.01925 

.09477 



.03136 

.05540 

4.01523 


.36081 

.23813 

. 13983 

.21390 


.01908 

.00025 

.01009 

.02335 



1.86352 


.41373 


.22663 

.14759 


.00779 

.17837 

.12503 




1.53699 

1.3096 

.05016 


.41466 

.26978 


.11014 

.16184 

.1187 


.00111 

.03102 

2.70158 
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TABLE B-6. U* Relative J'a Coefficients 

\ Final State energies are in cn . 

\ State 


Inltia^\ 

0 

289 

915 

1749 

2259 

44 21 

4585* 

5330* 

5596* • 

5791 

6283 

6445 

8394 

J8200 

.00650 

.00162 

.00192 


.00251 







.00056 


20SS1 

.1448 


.00981 

.00654 



.01436 

.00060 

.00155 



.00016 


21S16 

,01239 

.01020 

.00981 

.00654 



.00109 

.00419 

.00832 

.00026 


.00048 


22101 

.00032 

.00038 



.00100 

.00077 



.00027 



.00036 


2229B 

.00S2S 

.03019 

.01866 

.09120 

.00092 

.01361 

.00944 

. 00359 

.01165 

.00403 


.00108 


22780 

.03109 

.05402 

.05304 

.01134 

.00869 

.00863 

.01905 


.02838 

.00044 


.07049 


2286S 


.00381 

.00316 






.01258 


.02648 

.00124 


23370 

.06665 

.09078 

.03232 

.00532 

.01401 

.00252 

.00055 

.00859 

.01774 

.01560 

.00459 

.01381 


23779 

.01498 

.03009 

.03070 

.02228 

,00912 

.00634 

,00211 

.00182 

.01333 

.00082 

.00950 

.00210 

.00420 

24153 

.18125 

.1090 

.00837 

.03572 

.05346 

.03801 

.00557 

.00357 

.00258 

.00265 


.00824 


34367 

.0821 

.00289 

.03541 


.00460 

.00862 

.01216 

.00329 

.00146 


.01273 

.00412 

.00111 

24818 

.37751 

.19230 

.01024 

.02855 

.00768 

.00491 

.00909 

.00411 

.00281 

.01045 

.00158 



25680 

. 20692 

. 19050 

.10337 

.33885 

.04622 

.07231 

.02;-i)s 

.04532 

.04521 

.00518 

.00619 



26091 

.42247 

1.04236 


.36479 

.00795 

.00152 

.02936 

.00132 

.00333 

.00674 

.01059 

.00868 

.00199 

26303 


.60488 


.08387 

.02257 

.00493 

.02909 


.00639 

.00296 

.00279 

.00055 

.00093 

27088 

.00017 


.01953 




.00746 

.00525 

.00217 


.00,363 

.00782 


17357 

.06633 

.04732 

.23262 

.13637 

.08640 

.01499 

.02546 

.00518 

.07676 

.01550 

.00769 


.0054$ 

28042 

.25781 

. 12754 

.7673 

.03134 

.50432 

.05695 

.13488 

.00069 

.00571 

.00968 

.00260 



23614 

.05940 

. 04696 

.16110 

• .10838 

.08785 

.02228 

.02972 


.00680 


.00420 


.00220 

29504 

.03219 

.13796 

.06951 

.11384 

.06338 

.04923 

.01704 

,00098 

.03453 

.00491 

.01927 

.00013 


30061 

.08496 

.16710 

.08112 

.20851 

.C09S0 

.11560 

.12199 

.08448 

.01235 

.02739 

.05239 

.00028 

.00631 

30291 

.36416 

.12304 

.02672 

.18904 

.00907 

.00006 

.21983 

.19264 

.00416 


.00760 

.00030 

.00053 

30627 

.00143 

.11508 

.04146 

.00500 

.06259 

.00229 

.00276 

.04642 

.00927 

.00918 

.12093 

.00855 

.00609 

50901 


.16207 

.04607 

.22067 

.02117 

.01768 

.16549 

.13521 

.00184 


.02764 

.00027 

.00355 

31152 

.00271 

.36205 


.05456 

,03768 

.02186 

.00533 


.01552 

.02583 

.00304 


.08266 

32373 





.63559 


.00228 


.46165 

.25181 

.0a089 


.00272 

33290 


.19822 


.17152 

. 12050 

.01552 

.02651 

.10371 

.02197 

.01039 




33749 

.01140 

.13205 


.01000 

..01223 


.00458 



.CC2S0 

.01115 


,09405 

34471 

.09716 

.50653 


.14006 

.C68S2 

.02172 

.02135 

.06905 

.00656 

.04182 

.01115 


.09405 

34797 

1.0906 

1.50674 

.23678 

.54012 

.07210 

.00025 

.03849 

.09477 



.03136 


.05540 

35738 


.72161 

.23813 

.27966 

.42781 


.01908 

.00025 

.01009 

.02335 




37049 


.41373 


.22663 

.14759 


.00779 

.35774 

.29006 





35309 

*1.3096 

.05016 


.41466 

.26978 


.11014 

.16184 

.1U7 


.com 


.03102 
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APPENDIX C. THREE STATE RATE EQUATIONS 

The formalism for solving the rate equations for a three state 
plasma, will be developed here because it is used in the simple examples 


used in the text. 

For a three state gas the rate equations can be written 
dN, 

ar “ ■'u "2 " *13 *3 * *C1 - ^*21 * *31 * * 11 ’ "^l ' 

dN^ ^ 

a# = *21 **1 " *23 *3 * *C2 - (*12 " *32 * * 12 ’ *2 ’ 

dN„ 

dt^ = R 3 I Nj + R 32 N 2 + ^23 ^ ^3 > 

where R is the rate that state j goes to state i (both collisionally and 
radiatively) , R(-^ is the recombination rate into state i and Rj^ is the 
ionization rate from state i. 

In the steady state problem and when the total number of atoms of the 
species is known (N^) Eq. (C-lb) be rewritten 

‘^2l‘^^o ' ’ ^2^ ■" ^^23 ^3 " ^^2 ^ *^32 ^ ^2^ *^2 = *^02 . 'CC-2) 

combining terms and letting 


^2 *^12 * *^23 *^ I 2 * *^21 * 


(C-3a) 


®2 " ^23 " ^21 


(C-3b) 



^21 No 


(C-3c) 


Eq. (C-2) can be rewritten as 
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~^2 ^2 * ®2 ^3 " ~^2 * 

Similarily Eq. (C-lc) can be rewritten 

B3 N2 - A3 N3 = -C3 . 

Therefore, the state densities are 

* ^ 2^3 

2 = A2A3 - B2B3 ’ 

3 - A2A3 - B2B3 ' 


(C-4) 


(C-5) 


(C-6a) 


(C-6b) 


" ^^2 + N3) • (C-6c) 

To show how these are used let state 3 be the 3.9605 eV state in 

neutral uranium and state 2 account for all upper state interactions 

with state 3. Also the interaction with the ground state will be that 

of the 1.6691 eV state. The rates used will be for a neutron flux of 

2 X 10^^ neutrons/ (cm^-sec) and a temperature of 8000°K. Hence, the 

17 3 

total neutral uranium density is 2.58 x 10 atoms/cm . The other 

rate constants will be found in Table (C-1) and give values of 2.56 

X 10^^ atoms/cm^ for 1.15 x 10^^ atoms/cm^ for N2 and 9.63 x 10^^ 

3 

atoms/cm for N3. 



Table C-l. Rate constants for Eq. (C4) and Eq. (cg for 

14 2 

flux of 2x10 neutrons/ (cm -sec) and a 

temperature of 8000®K. 
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